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Abstract—The Terahertz (THz) band (0.3-3.0 THz), spans a 
great portion of the Radio Frequency (RF) spectrum that is 
mostly unoccupied and unregulated. It is a potential candi- 
date for application in Sixth-Generation (6G) wireless networks, 
as it has the capabilities of satisfying the high data rate and 
capacity requirements of future wireless communication systems. 
Profound knowledge of the propagation channel is crucial in 
communication systems design, which nonetheless is still at its 
infancy, as channel modeling at THz frequencies has been mostly 
limited to characterizing fixed Point-to-Point (PtP) scenarios up 
to 300 GHz. Provided the technology matures enough and models 
adapt to the distinctive characteristics of the THz wave, future 
wireless communication systems will enable a plethora of new 
use cases and applications to be realized, in addition to deliv- 
ering higher spectral efficiencies that would ultimately enhance 
the Quality-of-Service (QoS) to the end user. In this paper, we 
provide an insight into THz channel propagation characteris- 
tics, measurement capabilities, and modeling techniques for 6G 
communication applications, along with guidelines and recom- 
mendations that will aid future characterization efforts in the 
THz band. We survey the most recent and important measure- 
ment campaigns and modeling efforts found in literature, based 
on the use cases and system requirements identified. Finally, 
we discuss the challenges and limitations of measurement and 
modeling at such high frequencies and contemplate the future 
research outlook toward realizing the 6G vision. 


Index Terms—6G, channel measurements, channel modeling, 
channel sounding, deterministic modeling, diffuse scattering, 
molecular attenuation, propagation, radio frequency (RF), sta- 
tistical modeling, stochastic modeling, terahertz (THz), wireless 
communications. 
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I. INTRODUCTION 


S THE demand for higher data rates is growing expo- 

nentially, data rates that exceed | Tera-bit per second 
(Tbps) could be required beyond 2030 [1]. Very high data 
rates are anticipated in indoor locations where lower costs can 
be expected for extending wired backbone networks compared 
to using optical cables, while ITU [2] reported that by the 
year 2030, wireless systems will be as significant as wired 
cellular systems. A recently published white paper by Cisco, 
estimated that average Fifth-Generation (5G) cellular speeds 
will reach 575 Mbps by 2023, resulting in a 13-fold increase in 
the speeds recorded in 2018 [3]. Additionally, Ericsson high- 
lighted the imminent need for capacity, as global mobile data 
traffic is estimated to grow by a factor of 4.5 to reach 226 
ExaBytes per month in 2026 [4]. This unprecedented need 
for higher speeds and larger bandwidths in sight of the cur- 
rently saturated spectrum below 300 GHz [5], has motivated 
research in parts of the Radio Frequency (RF) spectrum that 
have been faintly explored, in order to meet the growing needs 
of the wireless communications industry [6]. A recently pub- 
lished white paper [7] derived within the “6G Flagship” by the 
university of Oulu in 2019, mentioned that Sixth-Generation 
(6G) systems should be able to deliver Tbps speeds and that 
the way to achieve these is through efficient utilization of the 
THz spectrum. THz will play a specially important role in 
6G networks, as it is envisioned that THz communications 
will empower real-life solutions and services in coexistence 
with the heavily saturated sub-6 GHz and millimeter-wave 
(mmWave) spectra, in order to meet key performance indica- 
tors such as peak rates of 1 Tbps, extremely Ultra Reliable 
Low Latency Communications (eURLLC) and high energy 
efficiencies of 1 Tb/J [8]. Henceforth, initiatives for realizing 
the 6G vision have been recently launched such as the “Next 
G Alliance” project by the USA and “Heat-X” European 6G 
flagship project [9]. 


A. Introducing the THz Band 


The frequency band that lies between the microwave and 
infrared regions of the electromagnetic spectrum spanning 
from 300 GHz to 3,000 GHz is termed as the THz band, with 
wavelengths ranging from | mm to 0.1 mm [10]. Recently, the 
THz band has received considerable attention as a promising 
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candidate for various applications that can be classified in the 
fields of wireless cognition [11], sensing [12], imaging [13], 
communications [14], and positioning [13], offering advan- 
tages such as ultra-fast data transmission [15], large band- 
width [16], and microsecond latencies [17]. THz propagation, 
has so far proven to be successful in spectroscopy for sensing 
applications such as material characterization [18], weather 
characterization [19] and study of astronomical objects [20], 
due to its frequency selective nature. Some first applica- 
tions identified where the THz wave could be applied to 
in 6G networks for wireless communications, include PtP 
scenarios such as kiosk downloads [21], wireless links in 
data centres [22], intra-device communications [23] and wire- 
less backhaul/fronthaul link extensions [24]. Nonetheless, as 
the technology continues to mature providing the appropriate 
beam-steering and tracking capabilities at THz, applications 
that include mobility can also be realized. Applications as 
such include Wireless Local Area Networks (WLANs), High 
Speed Train (HST) communications [25], and Vehicle-to- 
Everything (V2X) communications [26]. Last but not least, 
THz has found great interest in future nano-networks within 
the medical industry for in-body via nano-machine molecular 
communications [27], [28]. 

The aforementioned opportunities are ones that arise due 
to the unique characteristics held by the THz wave in addi- 
tion to the largely unoccupied and unregulated region that 
this end of the spectrum offers. For the development of an 
efficient communications system with high data rates, deriv- 
ing the characteristics of the propagation channel by accurate 
modeling is a necessity. So far, most studies have been limited 
to frequencies up to the mmWave spectrum (30-300 GHz) due 
to the lack of high frequency generators and sensitive detec- 
tors. At frequencies > 300 GHz propagation models need to 
consider high path loss, molecular absorption and diffuse par- 
ticle scattering (due to the short wavelength < 1 mm), as the 
most important phenomena that are very prominent, especially 
as we shift higher up in the THz frequency spectrum. 


B. Enabling Technologies for THz Waves 


THz research cannot be pondered as something “new”. The 
first paper ever published [29] dates back to 1893 and it 
involved THz spectroscopy. However, the term “THz” was 
not common and instead it was referred to as “obscure” radi- 
ation. In the experiments conducted there was a source and 
detector of THz radiation and by placing different materials 
in between, absorption coefficients were measured, much like 
a lot of the measurements conducted nowadays [30], [31], [32]. 
This shows that while a lot of work has been carried out and 
a lot of materials have been characterized, there is still a lot 
of interest and work to be done in basic spectroscopy at THz. 

What has happened recently, is that there is a growing 
realization that the photonics and electronics research coin- 
cides, and there is an ongoing merging of the electronics 
and photonics fuse of the world [33]. In recent years there 
has been a lot of progress in some key metrics, such as the 
power produced by various different sources vs frequency, 
where, e.g., the power produced by silicon sources in 2018 


SERGHIOU et al.: TERAHERTZ CHANNEL PROPAGATION PHENOMENA, MEASUREMENT TECHNIQUES AND MODELING 


has increased by many orders of magnitude compared to the 
power levels produced by silicon in 2008, while other exam- 
ples also show a similar sort of improvement [33]. This also 
shows that the production of integrated systems is possible 
and that enough power can be generated to do something 
useful with it. Various techniques are currently being carried 
out for the development of efficient sources and detectors of 
THz radiation [34], [35], [36], to bridge the ‘gap’ between 
conventional electronics and optics. Some of which, have 
made THz systems affordable and realizable by overcoming 
existing limitations through novel architectural designs and 
device structures. On the electronics side, solid-state elec- 
tronics technologies such as GaAs-based Schottky diodes, 
Monolithic Microwave Integrated Circuit (MMIC), InP het- 
erojunction bipolar transistors (HBTs), high electron mobility 
transistors (HEMTs) and Resonant-Tunnelling Diode (RTD) 
can now enable high power level generation from 100 uW 
up to the mW range [38], [37]. This shows an increase 
of more than double in a decade, especially at frequencies 
exceeding 1 THz. Likewise, various photonic technologies 
within the THz band have been developed recently, that 
enable W Continuous-Wave (CW) power generation, includ- 
ing photomixing approaches such as Uni-Travelling-Carrier 
Photo-Diode (UTC-PD) and also Quantum Cascade Lasers 
(QCLs) which are diode lasers capable of operating below 
2 THz [39], [40]. DFG/parametric generation within the 1-2 
THz range is also another option that enables mW radiated 
power levels to be obtained [41]. Technological growth in 
terms of power generation and detection therefore, seems to 
converge within the 3 THz expected limit. Thus realistically, 
work has to be done in that domain. 

On the other hand, there are many other important metrics 
besides power that require consideration, such as efficiency 
and sensitivity of detectors. These metrics are important as 
it is not just a matter of whether it is possible to gener- 
ate a THz signal but to also manage to tune and control 
things such as the direction, spectral content and polariza- 
tion of the THz beam [42], [43]. In the spectrum above 1 
THz, photonic based systems such as QCLs, Si-photonics and 
UTC-PD devices are expected to play an important part, espe- 
cially where tuning of > 1 THz is required [33]. It is equally 
important however, whether these things can be done fast and 
in a programmable and reconfigurable manner, which would 
require the use of phased arrays [44], [45] that are currently not 
commercially available. Phased array technologies that have 
recently emerged, include antenna arrays with a large number 
of elements, that can be placed compactly in a small amount of 
space to form Ultra-Massive Multiple-Input-Multiple-Output 
(UM-MIMO), based on an array with an ultra large number 
of elements, that can offer beamforming gains to overcome the 
severe path losses [46]. Reflective Intelligent Surface (RIS), is 
another enabling technology that differs from phased arrays, 
as it has no RF chain requirements, and allows optimization of 
the direction of the propagating wave that is directed upon it, 
by varying the Electro-Magnetic (EM) properties of the anten- 
nas via electronic control [47]. An RIS essentially would be 
employed to enhance the received signal when the direct path 
is not available. 


1959 


Metasurface architectures that form the basis for RIS, can 
also be employed to generate Orbital Angular Momentum 
(OAM) beams by exploiting the frequency selective nature 
of the THz wave. In particular, OAM can generate a great 
number of modes that are orthogonal to each other and 
which can enable Line-of-Sight (LoS) links diverted from 
spatial multiplexing, and enable OAM multiplexing in order 
to improve the throughput as an alternative to Orthogonal 
Frequency Division Multiplexing (OFDM), by alleviating 
interference issues to enable multiple access to users more effi- 
ciently, without consuming the available resources. In addition, 
spectrum access techniques such as Non-Orthogonal Multiple 
Access (NOMA) are also key THz enablers. NOMA exploits 
the fact that user experience will vary significantly due to 
the characteristics of the propagation channel that is affected 
by factors such as the mobility of users, frequency selec- 
tive fading and blockage effects, in order to provide more 
fair user access in terms of resource allocation, to enhance 
the overall networks’ spectral efficiency. However, the com- 
bination of NOMA with UM-MIMO or OAM, may only be 
desirable in LoS scenarios where spatial processing to sepa- 
rate the users doesn’t work and some gains can be achieved. 
Alternatively, the concept of Quasi-Optical Multiple-Input 
Multiple-Output (QO-MIMO) introduced in [48], overcomes 
interference issues associated with traditional Multiple-Input 
Multiple-Output (MIMO) systems by considering single-beam 
per feed antennas to exploit the multi-mode capacity of near 
field links. This method has shown to outperform traditional 
MIMO with large capacity gains and interference free signals. 
The disadvantage of approaches that aim to serve multiple 
users at THz, is that the aperture must be large enough com- 
pared to the Transmitter (Tx)-Receiver (Rx) separation for 
the channel to have enough ranks to support transmission 
of multiple parallel streams. Furthermore, the aforementioned 
technology enablers, introduce a number of resource man- 
agement challenges [8] that must be addressed individually 
and collectively to achieve the performance goals set. The 
performance must then be evaluated by an appropriate frame- 
work for RF/THz coexistence as well as THz-only operation in 
various configurations from which coverage probabilities can 
be accurately derived, to evaluate the performance of a user in 
a coexisting network, in order to maximize the coverage [49]. 
The framework developed must be aware of mobility chal- 
lenges faced in this hybrid deployment, in direct relation to the 
high molecular absorption and highly directive links expected 
at THz frequencies, that result in increasing its mathematical 
complexity massively [50]. Finally, Machine Learning (ML) 
can be considered as another key 6G enabler, where algorithms 
can be trained to offer intelligent mobility management solu- 
tions and efficient multiple access resource allocation, while 
also enhancing the diverse coexistence between THz and the 
lower frequency bands to guarantee the performance indicators 
promised [51]. 


C. Standardisation Efforts 


A first standard on wireless communications at 300 GHz 
has been defined by the IEEE 802.15 THz Interest Group 
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Fig. 1. 


IEFE 802.15.3d 252-321 GHz bandwidth allocation plan [54]. 


known as IEEE 802.15.3d-2017 [52], as an amendment to 
the IEEE 802.15.3-2016 [53] 60 GHz standard, for opera- 
tion up to 100 Giga-bits per second (Gbps). The 300 GHz 
system can utilize 8 different channel bandwidths between 
252-321 GHz at multiples of 2.16 GHz [54], as illustrated 
in Fig. 1. Both single carrier and on-off modulation schemes 
are being employed, which are combined with Low-Density 
Parity-Check (LDPC) and Reed-Solomon (RS) codes. Some 
of the applications identified for this range include PtP sce- 
narios such as kiosk downloads, data centre communications, 
fronthaul/backhaul, intra-device communications and WLANs. 

Spectrum allocation currently ends at 275 GHz for active 
services, but is not completely unregulated, as it is assigned for 
passive services such as Earth Exploration Satellite Services 
(EESS) and radio astronomy, that gather information based 
on the resonant frequencies of the molecules, and therefore 
these systems are only receiving and not transmitting [55]. 
These resonant frequencies are occupied within bands that 
cover most of the spectrum between 275-1,000 GHz, meaning 
that if we want to allocate these bands for active services and 
employ THz communications, this needs to happen by shar- 
ing the spectrum [56]. Furthermore, all frequencies between 
1,000-3,000 GHz may also be used for both passive and active 
services [56], though they have not yet been allocated to any 
service. In regard to the bands between 275-450 GHz, the 
ITU-R has invited its members to conduct sharing studies to 
see under which conditions and within which frequency bands 
sharing is possible. At WRC 2019 [55], four contiguous bands 
were identified between 275 and 450 GHz (275-296 GHz, 
306-313 GHz, 318-333 GHz, and 356-450 GHz) with large 
bandwidths, i.e., one of them with 94 GHz available for unre- 
stricted use in total, based on ITU-R studies [57], for land 
mobile and fixed services without the need to protect EESS. 
The Horizon 2020-EU-Japan joint project called ThoR, has 
also conducted a sharing study [58] for interference miti- 
gation, taking into account the operational characteristic of 
both passive and active services, and its findings are in line 
with the ITU-R studies. Three bands have been identified 
(296-306 GHz, 313-318 GHz and 333-356 GHz) in line with 
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bands where proper mitigation is required to avoid harmful 
interference and allow for coexistence with other applications, 
while for radio astronomy, on the contrary, proper mitigation 
techniques are already available [59]. 

In conclusion, Frequency bands > 275 GHz are inherently 
capable for implementing wireless systems with data rates 
of 100s of Gbps. Feasibility of wireless THz communication 
based on current developments and hardware demonstrations 
has been proven, while the first standard is available and so is 
the spectrum. Therefore, THz communication is a real option 
to be employed in future 6G wireless networks. 


II. MOTIVATION AND CONTRIBUTIONS 


In this section, the necessities revolving around channel 
propagation modeling and measurements at THz are high- 
lighted, and the existing works found in literature are com- 
pared and evaluated in terms of their contributions as well 
as their shortcomings. The research gap identified based on 
these shortcomings, is supported by arguments that provide 
work motivation and backing of the contributions offered in 
this work. 


A. Motivation and Related Works 


AS a prerequisite to communications system design and 
deployment, obtaining the characteristics of the wireless prop- 
agation channel between the Tx and Rx ends of the communi- 
cations system is imperative [60]. Knowledge of the properties 
and distinct features of the propagation channel in a given 
environment based on the frequency range of operation, allows 
performance limits to be set, that are crucial in the devel- 
opment of successful communications systems in terms of 
transceiver design, software architecture and protocols [61]. 
Therefore, studying the wireless propagation channel to ana- 
lyze its behaviour is critical in 6G communications system 
design. Obtaining the characteristics of the wireless propa- 
gation channel, necessitates correct procedures and empirical 
measurements to be conducted using channel sounders. The 
measurement data collected are then analyzed and modeled 
depending on the application scenario, and frequency range 
of operation. The applicability and performance of the models 
developed can then be compared as to which ones best describe 
the nature of the propagation channel in wireless communi- 
cations networks, while the characteristics obtained must be 
analyzed comparatively to enhance system design, and pro- 
vide a reference for future characterization efforts in higher 
frequency bands. 

Very interesting survey papers that cover channel mea- 
surements and modeling can be found in literature, that 
mainly focus on the 5G mmWave bands < 100 GHz. 
For instance, Wang ef al. [62] listed the requirements for 
modeling the mmWave bands, while the available mea- 
surement campaigns and models were reviewed. Similarly, 
Hemadeh ef al. [63], surveyed channel measurements and 
modeling efforts for the 5G bands, while challenges for 
system design and deployment were considered. Furthermore, 
Rappaport et al. [64] summarized measurements they have 
conducted at mmWave frequencies, at a period between 2011 
and 2013, while Salous et al. [65], presented an interesting 
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TABLE I 
TERAHERTZ PROPAGATION CHANNEL MEASUREMENT & MODELING WORKS IN THE LITERATURE 

Paper Title Year Source Type Paper Contributions Ref. 

1 Propagation Modeling for Wireless Com- | 2018 Article Channel models and considerations for SISO and MIMO [75] 
munications in the Terahertz Band systems operating in the THz band. 

2 Channel Modeling and System Concepts | 2020 Article Guidelines and considerations for fixed PtP scenarios in [76] 
for Future Terahertz Communications: the 300 GHz band based on the IEEE-802.15.3d standard. 
Getting Ready for Advances Beyond 5G 

3 Towards 6G: Paradigm of realistic tera- | 2020 Case Study A case study for realizing THz channels for close- [77] 
hertz channel modeling proximity communications and smart rail mobility. 

4 Survey on Terahertz Nanocommunication | 2021 Survey Review of applications, channel models in nano-networks 
and Networking: A Top-Down Perspec- and their requirements; discuss contributions toward dif- | [71] 
tive ferent layers of the protocol stack. 

5 A Tutorial on Mathematical Modeling of | 2022 Tutorial A tutorial on high-level mathematical framework for 
5G/6G Millimeter Wave and Terahertz 5G/6G including system level evaluation, mmWave and | [69] 
Cellular Systems THz KPI evaluation; Service processing at the BS side. 

6 Terahertz Wireless Channels: A Holistic | 2022 Survey A General overview of THz channel models, measure- [74] 
Survey on Measurement, Modeling, and ments and analysis up to 1 THz. 
Analysis 

This work | Terahertz Channel Propagation Phenom- | 2022 Survey/Tutorial| A survey/tutorial on propagation modeling at THz witha | - 

ena, Measurement Techniques and Mod- focus on 6G applications, particularly, using new absorp- 
eling for 6G Wireless Communication tion loss models, understanding channel sounding and 
Applications: a Survey, Open Challenges modeling for 6G applications, pointing out the need for 
and Future Research Directions AI to overcome the challenges faced. 


review on propagation phenomena and the progress of chan- 
nel sounding and modeling for mmWave channels. However, 
propagation characteristics differ from those found at THz due 
to the shorter wavelength and therefore so do measurement 
capabilities and modeling techniques. 

Many informative papers can be found that cover differ- 
ent aspects of THz communications, but which nonetheless 
discuss the propagation aspect briefly in relation to the envi- 
sioned applications. For example, Elayan et al. [66] reviewed 
the current technological progress at THz as well as general 
propagation characteristics and models, while Chen et al. [17] 
presented the technological progress along with envisioned 
applications at THz. Dang et al. [67] discussed potential appli- 
cation scenarios and services that can be provided for THz 
alone and also hybrid network operations. Tataria et al. [68] 
provided a discussion on the use cases and enabling tech- 
nologies that will form the basis of 6G networks such as 
new waveform design, UM-MIMO, novel signal process- 
ing techniques and core network design, while a tutorial 
presented in [69] describes a high-level mathematical frame- 
work for 5G/6G system design. In [70], the authors pro- 
vided a discussion on the defining features of THz wireless 
systems and four potential use cases related to 6G com- 
munications and the Internet of everything (IoE). Moreover, 
an overview for communications within nano-networks in 
the THz band was presented in [71], that covers different 
aspects of nano-communications and includes a summary 
of the available methods for modeling applications such as 
nano-device in-body communications and on-chip communi- 
cations. Some recent articles discuss channel measurements 
and modeling for 6G applications [72], [73], though they 
don’t focus on THz propagation characteristics, that vary sig- 
nificantly depending on the frequency range of operation, 
and peculiarities of the propagation environment. Therefore, 
there is a need to focus on the propagation aspect at THz 
frequencies, to demonstrate the applicability of the envisioned 
applications. 


A number of existing papers study THz channel prop- 
agation, providing general information on characteristics, 
modeling and measurements in the THz band. For example, 
Han et al. [74], presented a holistic survey on measurements, 
modeling and analysis in the THz band, however they focus 
on general characteristics of the propagation channel at THz, 
not analyzing comparatively the characteristics in relation to 
the envisioned applications, while they do not explore the 
potential the envisioned applications have in the atmospheric 
windows above | THz. An article presented in [75], consid- 
ers modeling methods at THz along with current issues and 
future considerations, though this is a general review paper 
that does not survey the literature comprehensively, so as to 
relate to potential 6G applications. Nonetheless, a few papers 
can be found that relate the propagation aspect at THz with the 
envisioned applications. A very useful article found in [76], 
provides guidelines and considerations for four potential PtP 
scenarios in the 300 GHz band, while a case study is provided 
in [77] for close-proximity and smart rail mobility scenarios; 
however they do not explore other potential applications, as 
well as the applicability and limitations of the applications 
considered above 300 GHz. 

Consequently, the existing survey articles do not ana- 
lyze comprehensively and comparatively the characteristics of 
the propagation channel, channel modeling and the require- 
ments for measuring at THz frequencies in relation to the 
applications identified, but cover different aspects of 6G com- 
munications at different frequency ranges where the channel 
characteristics and measurement capabilities vary significantly. 
In addition, they do not survey the literature comprehensively 
enough to discuss the feasibility of the envisioned scenarios 
in the THz regime, since they do not analyze the frequency 
bands above | THz and up to 3 THz where the technolog- 
ical advancements seem to converge, in order to investigate 
the feasibility of employing communications systems toward 
that end. Propagation related THz works are summarised in 
Table I, indicating a clear interest shift in recent years toward 
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the next generation of mobile communications. Much to the 
information provided on THz in light of the existing survey 
articles that cover different aspects of THz communications; 
an in-depth survey/tutorial in view of the current progress and 
future perspective of measurement and modeling of the THz 
channel at frequencies > 300 GHz for potential use case appli- 
cability in 6G networks, is still not available. In this survey, we 
aim to narrow this gap by addressing the shortcomings iden- 
tified and surveying the existing literature comprehensively 
and comparatively, aiming to provide a backbone that would 
motivate and help the research community in characterising 
the propagation channel effectively, for successful develop- 
ment and deployment of future 6G wireless communications 
systems operating in the entirety of the THz regime. 


B. Contributions of This Survey 


The aim of this paper, is to provide a thorough review of the 
propagation channel at THz for specific use cases and scenar- 
ios, giving recommendations that would enhance measurement 
procedures, and aid in the development of models that can 
characterize the propagation environment for the deployment 
of future 6G communication systems operating at frequencies 
between 0.3-3 THz. The main contributions of this survey are 
summarized as follows: 

e Potential applications that could be supported in 6G 
networks, and the requirements these applications pose 
on the underlying propagation channel in terms of mea- 
surement at THz, are discussed. 

e Propagation phenomena uniquely experienced at THz and 
their associated models are presented; Importantly, the 
absorption loss is computed up to 3 THz and total path 
losses arising from possible bands of operation for each 
application are assessed for the first time for 1-3 THz. 

e State-of-the-art channel sounding options at THz are 
reviewed, and compared for sounding the channel effec- 
tively in relation to the measurement requirements iden- 
tified for the 6G scenarios. 

e Stochastic, deterministic and combined stochastic- 
deterministic modeling approaches and their applicability 
at THz are assessed, for characterizing the propaga- 
tion channel accurately in relation to the identified 
applications. 

e Channel conditions resulting from the peculiarities of the 
propagation environment in relation to 6G scenarios are 
presented and analyzed. 

e Open areas of research are highlighted; In particular, 
pointing out the need to incorporate AI approaches in 
channel models to reduce the extent of measurement 
required, by capturing variations in space, time and 
frequency. 


C. Outline 


The structure of this paper is concisely illustrated in Fig. 2. 
The remainder of this paper is structured as follows; Initially, 
the use cases for wireless communications in the THz band 
as well as the system requirements for performing meaning- 
ful channel measurements are presented in Section II. The 
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Fig. 2. Structure of this paper. 


main propagation phenomena at THz and their modeling is 
presented in Section IV. THz channel sounding methods are 
introduced in Section V, which are necessary for conducting 
measurements based on the use cases identified in order to 
determine which method is best to be employed depending 
on each use case’s requirements. Channel modeling methods 
that can be followed for extracting the propagation channel’s 
characteristics in the THz band categorized as deterministic, 
stochastic and combined stochastic-deterministic are reviewed 
in Section VI. In Section VII, the most recent and signifi- 
cant efforts for characterising the propagation channel for 6G 
applications are summarised. Open issues and challenges being 
faced along with solutions and research directions for future 
modeling works are discussed in Section VIII. Finally, the 
findings of this paper are concluded in Section IX. 


III. USE CASES AND REQUIREMENTS 


For future communication systems deployment, the use 
cases identified lay a set of requirements and limits for 
performing meaningful channel measurement campaigns. On 
that account, this section presents the applications envi- 
sioned for THz communications whose applicability has 
been demonstrated in literature. Onward, a general set of 


TABLE II 
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TARGETED APPLICATION SCENARIOS FOR THZ COMMUNICATIONS AND MINIMUM REQUIREMENTS FOR CONDUCTING CHANNEL MEASUREMENTS 


Application Propagation Typical Propagation Character- | Antenna Alignment Re- | Ref. 
Environment Range istics quirements 
Kiosk Downloads Indoor; limited movement | ~ 1 m LoS, NLoS Manual alignment of Rx [143] 
by users by user 
Data-Center Networks Communications between | < 20 m LoS; D-NLoS Manual alignment [142] 
servers in data center 
Intra-Device Communica- | Communication within ~ 30 cm LoS; NLoS with signif- | Fixed alignment, low di- [23] 
tion devices; fixed links icant multipath contribu- | rectivity possible 
tion 
Backhaul/Fronthaul Fixed wireless link exten- | ~ 120 m LoS; Atmospheric attenu- | Highly directive antennas; [78] 
sions ation; fog and Rain atten- | fixed alignment 
uation; very high path loss 
WLAN Indoor; mobile users ~10m LoS; NLOoS; scattering Dynamic beam steering [79] 
HST Communication Communications within | <30m LoS; D-NLoS; | Dynamic beam steering; [154] 
trains and between station dynamically varying | high AoA estimation ac- 
infrastructure and train conditions; atmospheric | curacy 
attenuation; high path loss 
V2X Communication Communications between | < 20 m LoS; reflected path; dy- | Dynamic beam steering; [158] 
vehicles namically varying condi- | high AoA estimation ac- 
tions; meteorological at- | curacy 
tenuation; path loss 
Nano-Network Communi- | Nano-device < 10 mm Tissue penetration loss; | Fixed alignment [80] 
cation communication within scattering; molecular 
nano-networks; in-vivo noise 


system requirements for conducting practical channel mea- 
surements are described, by giving examples with regards to 
the 300 GHz band, in order to help in the selection of the 
appropriate measurement equipment for each of the use cases 
identified. 


A. Applications for THz Communications 


The governing factor that drives a meaningful measure- 
ment campaign, is the application scenario under investigation. 
Generally, all applications identified would benefit from THz 
links, that would enable high data rate transmissions. In 
addition, the spectrum above 300 GHz gives more licensing 
opportunities and flexibility than the lower frequency bands. 
Initially, it would be sensible to investigate scenarios that can 
be characterized by fixed PtP links where the Tx and Rx anten- 
nas are considered to be stationary for the whole duration 
of the measurement, such as intra-device, backhaul/fronthaul, 
data centre communications and kiosk downloads (though in 
this scenario the Rx side may be subject to minor displace- 
ment as the user would have to align the device to the Tx 
side). 

In the foreseeable future, advancing from fixed PtP scenar- 
ios to mobile scenarios, other applications are also envisioned 
in the THz band, such as communications within Vehicular 
networks, WLANs, and High Speed Train (HST) networks. 
These applications would benefit from MIMO capabilities to 
achieve extremely high data rates and also large capacities to 
accommodate a large number of users. Beam steering algo- 
rithms for device discovery requiring a high degree of Angle 
of Arrival (AoA) and Angle of Departure (AoD) estimation as 
the Tx and Rx will have no prior knowledge of each other’s 
position, would become essential in order to realize modeling 
of such applications. 


Table II establishes the requirements for channel modeling 
as to what exactly they are in terms of use cases, in relation to 
expected ranges and environments as well as expected propa- 
gation conditions and required antenna characteristics. The use 
cases identified for operation in the THz band are described 
as follows. 

1) Short-Range Communications: close-proximity commu- 
nications, include Kiosk downloads and on-board links, which 
may be considered as scenarios for file exchange between user 
devices and kiosks found in, e.g., airports, train stations as 
well as between devices such as routers and laptops on board 
a desktop. Propagation in such use cases may be assumed 
to be LoS over a short distance (< 1 m), with the poten- 
tial to provide “information showers” and achieve extremely 
high data rates targeting 100 Gbps [81] at 300 GHz. Thus, this 
prospective use case can be described as being highly applica- 
ble at frequencies > 300 GHz with the only uncertainty of user 
equipment antenna misalignment, where the user would have 
to essentially point the receiving device toward the transmitting 
end. 

2) Data Center Networks: Data centers are centralized 
units of operation, that house racks with servers that com- 
municate with each other by wired links [82]. A potential 
application of the THz wave at frequencies > 300 GHz 
because of the high bandwidth availability, flexibility it offers 
in hardware configuration/reconfiguration, lower costs as well 
as easier to maintain environmental conditions provided by 
cooling technologies for equipment operation. In addition, 
arrangement of racks in data centres is usually symmetric, 
and therefore propagation in such case may be LoS as well 
as Directed-Non-Line-of-Sight (D-NLoS) when the LoS is not 
available, over reasonably short ranges (< 20 m). 

3) Intra-Device Communications: The continuous rise 
in data transfer within modern day devices where 
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communications between components, different layers 
of motherboards and within on-chip networks are expected, 
has lead to complicated layout structures, reduced processor 
space and high energy level consumption [83]. As such, 
high data rate intra-device THz wireless communication 
links may prove to be a beneficial solution. The propagation 
environment to be considered for channel modeling is very 
small, in the order of a few cm due to the casing of the 
device. The antenna requirements may be more flexible with 
less directive characteristics but communications are expected 
to rely not only on LoS propagation but also by significant 
multipath scattering contributions within the small enclosure 
by Non-Line-of-Sight (NLoS) propagation. In this context, 
opportunities for this application are highly motivational for 
use in the bands higher up in the THz spectrum. 

4) Wireless Backhaul/Fronthaul: Wireless backhaul links 
are potential use cases for future THz communications, with 
LoS as the main propagation characteristic. A backhaul link 
is a fibre based or wireless link that connects the base station 
to the core network and is in the km range. Current wireless 
backhaul links are not sufficient for future capacity and latency 
requirements and therefore a potential solution might be that 
of THz wireless transmission as a wireless extension of the 
optical fiber between rural and highly populated urban areas 
to ensure high-speed connectivity everywhere in a cost effec- 
tive manner [84]. Current wireless backhaul implementations 
achieve up to 10 Gbps and the ThoR project [85] aims at real- 
izing data rates up to 100 Gbps for 300 GHz. THz Fronthaul 
links can be used to give broadband/high speed access to the 
end user in coexistence with the lower frequency bands in 
order to maximise the coverage, where the propagation is at 
LoS, and the link distance is much smaller compared to the 
backhaul. 

5) WLAN: Wireless connectivity that would require con- 
nections between various access points to replace fibre-optic 
cables and enable higher data rates. Typical indoor commu- 
nications in the THz band will involve scenarios such as 
office rooms and conference rooms. Such scenarios would 
have to consider LoS, link blockage, and scattering as the main 
propagation mechanisms and will be mainly applied for inter- 
connecting access points and user access at fixed locations. 
Communication links in such cases are limited by the enclosed 
propagation environment, with LoS and D-NLoS propagation 
mechanisms as well as scattering due to the enclosure, ranging 
up to a few meters. 

6) HST Communications: High data rate communications 
between trains, passengers and infrastructure to deliver high 
speed broadband access to HST users is expected to be deliv- 
ered by the next generation of mobile communications [25]. 
The most promising application of THz however, is in the 
Train-to-Infrastructure (T2I) channel to deliver very high data 
rates of 100 Gbps etc. Another potential scenario, is commu- 
nication within wagons for links between users and access 
points, where the large bandwidth at THz could be poten- 
tially exploited in the future for high data rate transmissions, 
and to accommodate a larger number of users in coexistence 
with the mmWave bands [25]. Propagation distances may vary 
depending on the scenario, ranging from a few meters within 
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the wagon up to a few tens of meters between infrastructure 
and train. Antenna requirements may also vary by scenario, 
where dynamic beam steering may be required to follow a 
moving train in the T2I scenario. 

7) V2X Communications: A network that is extremely 
intelligent and able to conquer the requirements of high 
data rates, ultra-fast, ultra-reliable and ultra-low latency 
exchange of information, is a future vision for 6G V2X 
communications [26]. The need for enhanced road safety 
and inter-connectivity that satisfies the requirements men- 
tioned, can be aided by communications in the THz band, 
to enhance the throughput. Particularly, highly directive THz 
links in the lower end of the spectrum may be utilized for 
large information exchange within vehicular fleets of semi- 
autonomous vehicles performing cooperative driving [87] as 
well as for data relaying in congested networks to enhance 
connection reliability [88]. At THz it may be necessary to 
have such communications between vehicles at points where 
they are stationary or moving slowly in a short range, such as 
to do ‘vehicle kiosk downloading/uploading’ with LoS propa- 
gation between two vehicles in line with NLoS contributions 
and scattering from nearby buildings and vehicles. 

8) Nano-Networks: The THz channel has been found to be 
a candidate for many applications including those in the med- 
ical domain. Nano-machine enabled communications such as 
in-vivo communications within sub-mm sized regions inside 
the human body that are otherwise not accessible by any other 
medical device. For modeling the channel for in-vivo appli- 
cations, the human tissue needs to be considered. Path loss 
can be computed based on different molecular characteristics, 
scattering from cells within the tissue, distance between the 
nano-machines, molecular noise and frequency of operation. 
Nano-network communications can be considered as highly 
applicable for operation in the entirety of the THz band, due 
to the short communication distance, and wavelength compa- 
rable to that of different molecules present in human tissue, 
that would allow for the identification of spectra that cannot 
be captured otherwise in the lower frequency bands. 


B. Requirements for THz Channel Measurements 


1) Bandwidth: While the spectrum is vastly available 
above 300 GHz in the THz regime, channels at multiples 
of 2.16 GHz have been introduced that allow for contigu- 
ous bandwidths of a few 10s of GHz targeting data rates of 
100 Gbps. More specifically, as per [52], data rates around 
100 Gbps can be achieved with an 8-PSK modulation scheme 
with about 50 GHz bandwidth. Higher order modulation 
schemes should be employed to give higher data rates, such 
as 64-QAM that is able to provide data rates > 200 Gbps for 
the same bandwidth (~ 50 GHz) [52]. To satisfy the capacity 
requirements and achieve high data rates, channels composed 
of many 10s of GHz of contiguous spectrum that can be freely 
allocated for active services with high order modulation and 
sophisticated coding schemes will be required. 

2) Antenna Directivity and Beamwidth: Side-lobe suppres- 
sion and spatial filtering are the main requirements for anten- 
nas operating at THz, something that cannot be provided 
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by omni-directional antennas used in the lower frequency 
bands. Therefore, highly directive antennas with high gain 
are required, which are very much necessary to combat the 
high path loss and large atmospheric absorption effects of 
the highly frequency selective THz channels and also pro- 
vide strong enough signal components that would be detected 
by the receivers. Antenna directivity can be described as the 
ratio of the maximum power density to its average value over 
all directions [89]. When no information is present on the 
effective aperture, the directivity in relation to the antenna 
beamwidth can be approximated by: 
41253 

Pup PHP 
where f;p is the half-power beam width in one principal plane 
and @7;p is the HPBW in the other principal plane, both mea- 
sured in degrees. HPBW is the angle in which the radiation 
intensity is at least half of the peak radiation intensity. For 
example, considering an HPBW of 10° on both planes, that 
would give us an antenna directivity from the radiation pattern 
of ~26 dB above isotropic. The gain of the antenna can then 
be found by multiplying the directivity value with the radia- 
tion efficiency of the antenna. In presence of multiple antennas 
that form a phased array and which are required for dynamic 
beam steering, planar arrays are the preferred solution at THz 
due to the advantages they hold in terms of large beam form- 
ing gains, fast beam-steering and side-lobe attenuation [90]. 
The directivity of a planar array including information on the 
beamwidth can be approximated as follows: 


32400 
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the importance here is to incorporate antenna patterns into 
channel simulators, so that they are accounted for in the 
capacity/data rate computation. Then patterns can be done by 
trigonometric functions to form the main beam, and an angular 
mask to cover the peak of the side lobes at any given angle. 

3) Link Distance: The measurement distance separating the 
Tx and Rx sides of the communications channel, determines 
the extent to which accurate channel measurements at THz 
frequencies can be conducted. Indoor measurements require 
Tx-Rx separation distances of a few meters, while in outdoor 
scenarios 10s and even 100s of meters of measurement dis- 
tances will be required, where various effects would have to be 
catered for such as a large number of reflections from build- 
ings, vegetation and fading from meteorologically induced 
attenuation. Cabling is a largely negative factor that currently 
affects measurements over long distances. A Radio Frequency 
Over Fiber (RFoF) extension for frequency channel sounding 
was developed for the first time in [91] to extend the Tx-Rx 
channel measurement distance above 100 m and was success- 
fully demonstrated in an outdoor environment. An important 
factor here is the system limitation when a THz radio link has 
several GHz of system bandwidth. The standard definition for 
noise figure is derived from the added noise by the receiver, 


NRec: 
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where No is the noise power added by the receiver, and 
both have corresponding noise temperatures from use of kTB. 
Therefore taking an ambient temperature To of 290 K, this 
calculates to form —84 dBm for 1 GHz bandwidth then up 
to —74 dBm for 10 GHz bandwidth. This noise floor would 
increase by up to another 10 dB based on the noise figure at 
the receiver, defined as: 


T, 
NoiseFigure = 10log19 (1 + Tre) (dB) (4) 
0 


where TRec is the temperature of the added noise from the 
receiver. So what is measured over fibre might exceed what is 
achievable actually by a THz system. 

4) Dynamic Range: The dynamic range of the equipment 
used for measurements is very important in determining the 
signal strength of different Multipath Components (MPCs). 
A high order of diffuse reflections coming from the very 
short wavelength of the THz wave interacting with irregu- 
larities on surfaces similar in magnitude to the wavelength, 
and indeed significant contributions of specular components 
resulting from the highly directive beam are expected. For 
that reason, the dynamic range of the measuring equipment 
will help in identifying these resulting components, that will 
be very helpful to determine the effect they have on the over- 
all received signal strength. Therefore, the noise floor must 
be at a low enough level to achieve a high dynamic range, 
quantified as the difference between the strongest power com- 
ponent received and noise level. For the SC channel sounder, 
the dynamic range can be calculated as a function of the 
m-sequence length calculated as: 


Dynamic Range = 20logj9(m) (dB) (5) 
where m, would have to be significantly large at frequencies > 
300 GHz, to aid in the detection of much weaker multipaths, 
especially for applications that suffer from high atmospheric 
absorption due to longer communication distance require- 
ments, in order to combat the problem of channel sparsity 
that is prominent at THz. Generally, all applications identified 
would require a high dynamic range; however, applications 
which may be less affected by the high losses present in 
the atmosphere such as intra-device links and short-range 
communications, a lower dynamic range should suffice. To 
achieve a good dynamic range of 80 dB, m must equal to 
10,000. Dynamic range can be further improved by reducing 
the noise floor through averaging or by IF bandwidth reduc- 
tion, in exchange for longer measurement times. Another, but 
more complicated option that can improve the dynamic range 
of channel sounding systems, involves employing the RFoF 
technique due to the lower losses incurred compared to RF 
cabling. 

5) Doppler Sampling Rate: The Doppler effect, is where 
the propagating waves’ frequency changes based on the veloc- 
ity of a moving object in relation to the sources vicinity. E.g., a 
user moving across the room and blocking the signal transmis- 
sion in an indoor scenario, or equivalently a vehicle moving 
away from a stationary source located on another vehicle in a 
nearby lane in an urban scenario. The need to measure Doppler 
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spectra, is down to the impact the Doppler has on the wave- 
forms used, such as OFDM with multiple carriers, where too 
much Doppler spread at THz level is not conducive, since it 
can cause loss of sub-carrier orthogonality and can limit the 
achievable data rates [92], [93]. Weakened sensing accuracy of 
multiple carrier systems, is also another limiting factor induced 
by high Doppler spreads, that also affects the data rates unfa- 
vorably [94]. To measure the Doppler spectra in a static case 
practically, we should rather evaluate the maximum Doppler 
shift in fractions of wavelength moved. Assuming a frequency 
of 300 GHz, this gives a wavelength of 1 mm. We can move a 
sliding sample by As, which corresponds to a change in spec- 
ular or diffuse reflection position Ap that is caused due to 
the roughness of the surface. We need to ensure Ap is small 
enough that it does not move beyond a portion of a wave- 
length from sample to sample. We can form a derivation from 
the Doppler shift as follows: 


Up Up CO Up 

k= Mash P habe NT (6) 
thus, for the representative velocity of the specular point, vp, 
this relates to how we sample relative to a wavelength, as 
we can say that the sampling frequency equals to fg = kfg. 
Therefore as we want to increase the sampling frequency rela- 
tive to the Doppler frequency by factor k, where 4 is a suitable 
choice in accordance with the Nyquist-Shannon sampling the- 
orem [95], then we need to sample every quarter wavelength. 
Thus, Ap must not move more than 0.25 mm at every sam- 
ple point. This depends on how rapidly Ap changes as As 
changes. With reasonably ‘smooth’ surfaces such as wood, we 
can likely expect 0.5 mm resolution is sufficient for coherent 
sampling in order to avoid aliasing effects, with a high auto- 
correlation above 0.9 being desirable, that must be checked 
following measurement. This means that measurements at THz 
frequencies may be limited by substantial practical difficul- 
ties requiring a trade off between accuracy and feasibility, or 
otherwise, requiring costly high precision positioning equip- 
ment. A similar process in measurement as such, has been 
recently adopted in [96] in order to observe PtP variations 
in multipath distribution by reflections off indoor surfaces. 
In the case where the Doppler spectra are measured using a 
sliding correlator dynamic channel sounder where we essen- 
tially transmit a maximum length, m-sequence, the maximum 
Doppler frequency can be obtained as 1/2km7 where k repre- 
sents the scaling factor of the sliding correlator, and 7 is the 
clock period, leading to, 


a Er (ms) o) 


meaning that the velocity of the observer is inversely propor- 
tional to m, suggesting that increasing its value can enable 
longer time delays to be resolved, while lowering the velocity 
would allow greater Doppler shift resolutions. For example, to 
measure the Doppler shift of a user blocking the LoS signal 
in an indoor scenario, given a standard value of k = 5000, 
0.1 ns and a user velocity of 7 km/h at 300 GHz, 
we would require an m-sequence length of 514. These fac- 
tors would depend on the capabilities of the dynamic channel 
sounder under consideration. 


T = 
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6) MultiPath Resolution: High delay resolutions that would 
enable multipaths to be resolved by their Time of Arrival 
(ToA) in the order of the wavelength are expected at THz, 
due to the large bandwidth availability, something which was 
not possible in the lower frequency bands. The delay res- 
olution can be calculated as 1/bandwidth. For instance, at a 
centre frequency of 300 GHz, this means we have a wave- 
length of 1 mm which corresponds to a difference in delay 
between multipaths of 3.33 ps. In order to capture channel 
variations and observe components arriving in the order of 
the wavelength, we would require a bandwidth of 300 GHz 
(150-450 GHz). For dynamic measurements, the clock rate of 
the channel sounding device determines the delay resolution 
offered by the system. This parameter is extremely important 
when performing channel measurements in environments that 
are heavily influenced by scatterers such as WLANs where 
a high delay resolutions would be required, in antithesis to 
urban environment for applications such as fronthaul/backhaul 
where LoS is expected to be the main propagation 
mechanism. 

7) Other Requirements: The choice of channel sounder to 
be employed, largely depends on the given application sce- 
nario e.g., whether its for material characterization, a fixed 
PtP communication scenario or a mobile scenario and what 
other requirements the targeted application has, such as band- 
width, speed, accuracy etc. What is a common need, is 
that of a relatively cheap and compact channel sounder that 
is easy to transport and one where cabling between the 
channel sounder and Tx/Rx would be as less complicated 
as possible to avoid calibration issues that are very com- 
mon when performing these measurements. Furthermore, for 
short range scenarios, the channel sounder itself when large 
in size may act as an undesired reflector. In addition, to 
allow for accurate THz channel measurements, the channel 
sounding system has to be characterized itself in order to estab- 
lish metrological traceability to the International System of 
units (SI) [97]. 


IV. THZ PROPAGATION PHENOMENA 


Accurate channel modeling and characterization is required 
for successful deployment of THz systems that can comply 
with industrial demands. To characterize the THz channel 
accurately and efficiently, it is of high importance that the 
nature of the propagation environment at this frequency range 
is considered. The THz wave possesses many unique charac- 
teristics that make it advantageous in comparison to microwave 
and mmWave communications, such as increased speed and 
capacity. Propagation phenomena considered in the microwave 
and mmWave bands include path loss, shadowing, penetration 
loss and diffraction. In the THz band, different propagation 
conditions besides the high path loss become much more pro- 
nounced such as molecular absorption and particle scattering. 
The effects of rain, fog and cloud also need to be taken into 
account as in the mmWave bands, while the effects of dust [98] 
and turbulence [99] are also limiting propagation phenomena, 
though they are often not as significant. 
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A. Beer-Lambert Propagation Model 


As the wavelength of the THz wave approaches the 
dimensions of molecules found in the atmosphere and in 
human tissue, this results in strong absorption and scattering 
of the wave. The High Resolution Transmission Molecular 
Absorption (HITRAN) spectroscopic database [100], is the 
most commonly used tool that provides a set of parameters 
that are employed in different modeling codes, to compute 
the attenuation experienced by the propagating wave. A useful 
series of tutorials provided by the creators of HITRAN on how 
to extract spectroscopic data of absorption begins in [101]. 
Computation of the absorption coefficient is based on the well 
known Beer-Lambert law of attenuation [102], that takes the 
given cross-sections of a gas that are found empirically by 
varying the temperature and type of molecule, and calculates 
the absorption based on the wavelength, gas concentration, and 
path length through the gas. The contributions of total molec- 
ular absorption are dependant on each isotopologue from each 
gas present in the medium. The absorption cross-sections Ky 
(cm?/molecule) can be calculated by: 


l 
pa (cm? /molecule) (8) 
pL 
where 7, is the spectral transmittance (cm7!) 


(molecule x cm7?) at a wavenumber v, temperature T (K) 
and pressure P (pa), where p is the density (molecule/cm*) 
along the path length L (cm). The absorption coefficient 
is then calculated as the absorption cross-section times the 
number of molecules in the path. The molecular specific 
attenuation can then be calculated as 20log(exp(a)), where 
a is the absorption coefficient. The modeling codes that are 
then used to compute various attenuations, require application 
of radiative transfer theory and input characteristics such as 
line shapes, scattering, continuum absorption etc. [103], that 
increase the complexity of the channel model and are gener- 
ally challenging for the beginners of the HITRAN database. 
The total path loss experienced by the system is then found 
as the sum of the specific attenuation due to absorption and 
Free Space Path Loss (FSPL) experienced in the medium it 
interacts with. An additional effect experienced, is that of 
molecular noise that causes molecular vibrations excited at 
the resonant frequencies of the molecules, causing emission 
of EM radiation, degrading the signal strength [104]. In terms 
of propagation in EM nano-networks, the total system noise 
takes additional contributions from the receiver antenna and 
from devices that are located in the immediate vicinity of 
the source. The concept of molecular absorption needs to 
be heavily accounted for, in networks operating under the 
NOMA principle. Novel User-pairing schemes and sophis- 
ticated channel models based on the complex mathematical 
structure of the Beer-Lambert law which must adapt to these 
characteristics according to the effect of molecular absorption, 
are therefore necessary [105]. 


B. Free Space Path Loss 


FSPL, is defined as inevitable attenuation experienced by 
a wave propagating in free space between two isotropic 
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Fig. 3. Free space path loss for the frequency range 1-3,000 GHz with unity 
gain Tx and Rx antennas. 


antennas [106]. It is described by the Friis’ transmission 
formula given by Eq. (9), 


PA? (4r)? d? Aet Aer 


= P; Act Aer 
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P,(d) = (9) 
where A is the wavelength, d is the distance the wave will 
travel, with P, being the transmit power and Aet, Aer being 
the Tx and Rx antenna effective areas. Based on Eq. (9), Fig. 3 
shows the FSPL for a PtP link for three different separation 
distances up to 3 THz. Initially, substantial path loss well 
over 100 dB is formed at 100 GHz, with a 20 dB per decade 
increase resulting in an increase by 20 dB from 0.1 THz to 
1 THz and nearly 10 dB increase from 1 THz to 3 THz. 
However, in Fig. 3 it is assumed that the antenna gains are 
equal to unity and the attenuation observed is simply the path 
loss. If the gain includes information on the aperture area of 
the Tx (Aet) and Rx (Aer) antennas, meaning that if the ratios: 


Act and Aer 
À À 
are maintained at any frequency, the received power P,(d) 
decreases with the square of the frequency and distance, 
i.e., the FSPL attenuation increasing proportionally with the 
decrease in P,(d), and hence, the path loss does not change. 
Therefore the antenna size requirements as we go up in 
frequency reduce. But at the same time it is generally not pos- 
sible to increase Aet or Aer at any single frequency beyond 
a corresponding limit, because the aperture efficiency falls 
consistently in cases such as horn antennas or parabolic dish 
reflectors. Array antennas as the other solution also have an 
equivalent limitation in that the increase in number of ele- 
ments will result in higher feed losses. Also at THz the closer 
physical spacing where feeds do not physically size down will 
cause substantial practical difficulties and further lossy charac- 
teristics. Finally, an omnidirectional antenna cannot increase 
its aperture, which is already nearly as low as an isotropic 
antenna, by virtue of its design. Therefore, based on the above 
observations, though in theory an increase of Aet and/or Aer 
will overcome path loss, it is inhibited by physical limitations. 


1968 


IEEE COMMUNICATIONS SURVEYS & TUTORIALS, VOL. 24, NO. 4, FOURTH QUARTER 2022 


100,000 i | 

| ll 

_ 10,000- : 

E EER een eee, OO Gd | OR a el Fy Y. 
+----f[- --ff----#-p---]-§ 44 - f--- My 

S 1,000+ sa <2 2000 dB/km 

T ; <> wii v2 

Cc 

S í W7 W8w9 W10 

S 1o0Ff----14W---Hw6--Ẹ}------4------4------4------4------ -1 - 

3 L <> 

2 N m Lae 100 dB/km 

E 1 3 

<x ! A 

2 1 EN) i E 

5 7. 

Fi ef 

Q I 

D 1 1 + 

£ 1 

= | 300 GHz 

0.1 | = Water vapor absorption (T= 15°C P=101.3 kPa p=7.5 g/m’) 7 
1 =—=Dry air (oxygen absorption) 
0.01 i | | | fi 
0 500 1000 1500 2000 2500 3000 


Frequency (GHz) 


Fig. 4. Atmospheric attenuation in free space for the range 1-3,000 GHz. 


C. Attenuation by Atmospheric gases 


At THz the main form of attenuation of the signal is 
expected due to meteorological attenuation known as atmo- 
spheric absorption. When the wavelength of the THz wave 
approaches the dimensions of water and oxygen molecules, 
resonances are generated at selective carrier frequencies result- 
ing in large attenuation peaks that degrade the signal severely. 
In Fig. 4 we can observe the atmospheric attenuation by 
water and oxygen molecules over the 1-3,000 GHz frequency 
range. It can be observed, that the contributions of oxygen 
for the bands above 300 GHz are insignificant and the atmo- 
spheric attenuation is dominated by water vapour. To calculate 
the atmospheric absorption at THz, it is necessary to extract 
spectroscopic data of absorption from the HITRAN database 
and apply radiative transfer theory. The computation of the 
atmospheric absorption can be done by first calculating the 
absorption lines of each gas molecule. When gas molecules 
collide with each other, the effect of Lorentz line broadening 
forms, that can be described by a half-width parameter that 
defines the width and shape of the lines, which can be com- 
puted as the average time between collisions, and which varies 
with temperature and pressure. The Van Vleck-Weisskofp 
(VVW) line shape [107] is employed to model the shape of 
the absorption lines. Then, based on the Beer-Lambert law, the 
molecular absorption spectrum is computed. Alternatively, the 
atmospheric attenuation for 10-1,000 GHz can be calculated 
by the ITU-R p.676-10 model for attenuation by atmospheric 
gases [108]. 

As illustrated in Fig. 4, below 300 GHz the attenuation is 
minimal reaching up to a few dB while in the underlined spec- 
tral windows below 1 THz, W1 = [331-375 GHz], W2 = [385- 
440 GHz], W3 = [453-510 GHz], W4 = [630-705 GHz], and 
WS = [815-890 GHz], the attenuation ranges up to 100 dB/km 
which still is a substantially massive amount of attenuation 
to achieve long distance communications, comparable with 
the attenuation we expect from a Faraday cage. In addition, 


different frequency-selective attenuation peaks are observed 
due to the resonant behavior of molecules in the atmosphere, 
causing severe attenuation levels some of which exceeding 
10,000 dB/km that would have to be strictly avoided. In 
the region between | to 3 THz the signal undergoes signif- 
icant attenuation by water molecules in the atmosphere and 
the next atmospheric windows with attenuations below 2000 
dB/km are found in, W6 = [995-1078 GHz], W7 = [1241- 
1320 GHz], W8 = [1325-1394 GHz], W9 = [1430- 
1589 GHz], W10 = [1815-1848 GHz], W11 = [1938- 
2013 GHz], W12 = [2088-2136 GHz] and W13 = [2495- 
2562 GHz] that may be potentially utilized only for near-field 
communications (< 0.1 m) [109]. 


D. Rain and Snow Attenuation 


Additional forms of meteorological attenuation that the THz 
wave may experience when propagating in an outdoor envi- 
ronment are those of rain and snow [110]. The ITU model 
provided in [111], describes the rain attenuation, yp in dB/km 
experienced by the THz wave over the 1-1,000 GHz frequency 
range given by the power-law relationship in Eq. (10), 


yp =kR* (dB/km) (10) 


calculated by the rain rate R (mm/h) where k and a are 
frequency dependent coefficients derived from scattering cal- 
culations at vertical and horizontal polarizations. As the rain 
droplets’ shape and size are in similar order to the wavelength, 
the propagation losses must be represented as a function of 
polarization and wavelength. In Fig. 5 the rain attenuation for 
horizontal polarization regarded as parallel to ground propa- 
gation for the 1-1,000 GHz frequency range can be observed, 
representing the worse case of attenuation. Four different sce- 
narios of rain are shown where the rain rate can vary from 
less than 1 mm/h for light rain to over 50 mm/h for extreme 
rain conditions. These losses are seen to flatten out around 
100 GHz, beyond which no additional attenuation is expected 
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Fig. 5. Rain attenuation for horizontal polarization in the frequency range 
1-1,000 GHz. 


within the bands of interest up to 3 THz [112]. The reason for 
the constant attenuation at higher frequencies is that particle 
scattering has moved from Rayleigh to Mie scattering [113] 
due to the shorter wavelength making it more applicable, indi- 
cating that particle scattering with pattern similar to that of 
a dipole radiation pattern (Rayleigh scattering) changes to a 
more direct scattering pattern with a forward lobe (Mie scatter- 
ing). These losses may not be as significant as the path losses 
and atmospheric losses but will have to be considered regard- 
less for accurate channel characterization. With regards to 
attenuation due to snow, an equivalent ITU propagation model 
does not exist and the snow particles are merely considered 
as rain droplets [114]. 


E. Cloud and Fog Attenuation 


Another meteorologically induced condition similar to that 
of rain and snow that can cause fading in an outdoor envi- 
ronment, is that of fog and cloud. Fog and cloud may be 
considered as the same atmospheric phenomena, that vary only 
by height above ground. Both phenomena are defined by the 
ITU model [115] under the assumption that the wave trav- 
els solely through a uniform fog or cloud environment. The 
specific attenuation ye can be computed by Eq. (11), 


(11) 


where K; is the cloud water specific attenuation coefficient 
((dB/km)/(g/ m*)) found by applying a Rayleigh based double- 
Debye scattering model that is used to extract the dielectric 
permittivity of water and which holds for frequencies between 
1-1,000 GHz. M (g/ m) is the liquid water density in the 
fog or cloud while T and f are the temperature and car- 
rier frequency, respectively. In Fig. 6 we can see the specific 
attenuation effects of fog/cloud at a temperature of 20°C for 
1-1000 GHz carrier frequencies that are typically around 0.05 
(g/ m*) for medium fog/cloud where the visibility range is in 
the order of 300 m, and 0.5 (g/ m?) for heavy fog/cloud where 
the visibility range drops to 50 m. The losses, similar to those 
of rain and snow, are much smaller than those of path loss and 


ye(f,T) = Ki(f,T)M (dB/km) 
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Fig. 6. Fog or cloud specific attenuation for the range 1-1,000 GHz. 


atmospheric absorption, and can be easily compensated for, as 
the attenuation forming at 100 GHz is around 5 dB/km with an 
increase of about 15 dB from 0.1 to 1 THz for the worst case 
condition of heavy fog/cloud. Beyond this region and up to 
3 THz, the attenuation formed does not vary significantly and 
increases by just a few dB [112] though no standard currently 
exists. 


F. Earth Attenuation 


The presence of Earth in an outdoor environment plays a 
significant role in signal attenuation at THz and may reduce 
the maximum achievable data rates [116]. The electrical char- 
acteristics of the surface of the earth, including pure water, 
sea water, ice, soil and vegetation cover, for frequencies up 
to 1,000 GHz can be computed using the ITU-R P.527-5 rec- 
ommendation [117]. High amounts of absorption would be 
expected due to the high water content, while the effects of 
scattering need consideration as to whether earth components 
can contribute to the received signal by NLOS propagation. 
Currently, no specific model exists that computes the propa- 
gation loss due to earth other than the ITU model in [118] 
for vegetation attenuation, which holds for frequencies up 
to 60 GHz. This is because the lack of measurement data 
makes it difficult to develop an appropriate model at THz 
frequencies. The THz-TDS technique could be utilized to mea- 
sure the water content on leaves and soil, and the data collected 
could be later on fed into models that characterize the outdoor 
propagation environment [119]. 


G. Propagation Mechanisms 


LoS can be considered as the most prominent propagation 
mechanism at THz as a result of the highly directive pen- 
cil beams expected to be transmitted by antennas operating 
at these frequencies. D-NLoS propagation, that was initially 
proposed in [120], is considered as a useful mechanism that 
would enable communication links when the direct LoS is 
not present, empowered by specular reflections off surfaces 
with good reflective properties as well as by RISs, to direct 
wave propagation in more complex cases over different angles. 
D-NLoS differs from the traditional concept of NLoS that is 
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TABLE III 
CATEGORIZING PROPAGATION LOSSES BASED ON ATMOSPHERIC WINDOWS AND EXPECTED LOS RANGES IDENTIFIED FOR EACH USE CASE 


ius Propagation ~ Expected ~ Propagation Losses (dB) 
Ref. Use Case/Application Phenomena Ranges 
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Ht 4 0 = = = B B È FB B FSF FE BF FB B 
[23] Intra-Device @®@@® O O 0.3 m 0.05) 0.05) 0.07| 0.07| 0.07| 0.11| 0.19| 0.23} 0.21} 0.2 | 0.33| 0.53| 0.57 
[143] | Kiosk/Desktop @®@o O 1m 0.16} 0.17) 0.23} 0.24] 0.24) 0.36) 0.65} 0.75} 0.69) 0.66} 1.09) 1.77] 1.91 
[79] WLAN e © O O 10 m 1.6 | 17 | 23 | 24| 24] 35 | 65] 7.5 | 69 | 66| 10.9) 17.7) 19.1 
[142] | Data Centre @®@eo°o 20 m 3.2 | 34] 46] 48 | 4.80) 7 13 | 15 | 13.8) 13.2) 21.8) 35.4] 38.2 
[158] | V2x © © © © 20 m 3.67| 3.96| 5.15| 5.46| 5.54) 7 13 | 15 | 13.8} 13.1] 21.8| 35.4| 38.1 
[154] | HST © © © © 30 m 5.5 | 5.9 | 7.7 | 8.2 | 8.3 | 10.6| 19.5) 22.6| 20.6| 19.7| 32.6| 53.1| 57.2 
[78] Backhaul/Fronthaul © © © © 120 m 22 23.8) 30.9| 32.8) 33.2) 42 78 90 82.8} 78.6} 130] 212| 228 


mainly used for sub-6 GHz communications as non-directed 
NLoS contributions to the received signal strength may be 
in many cases insignificant at THz due to high atmospheric 
losses, as they would be lost in thin air. However, for cases 
where the propagation environment is confined to a small 
space, such as in the case of intra-device communications, 
the NLoS components are expected to play a major role. 

Similar to both microwave and mmWave radiation, THz 
waves can penetrate through different non-conducting materi- 
als such as plastic, clothing and paper [121]. However, strong 
absorption of the wave, depending on material structure and 
thickness is extremely possible. Penetration at THz is lim- 
ited by the low output power levels and varies by thickness 
and composition of the material it interacts with [122], The 
THz wave is non-invasive to the human body, and is therefore 
highly susceptible to link blockage in indoor locations (i.e., 
when a user is moving across the room) which becomes a 
major issue when employing communications systems with 
highly directive antennas, hence, the need for propagation 
mechanisms such as D-NLoS. Furthermore, diffraction effects 
common in the lower bands < 100 GHz that can be described, 
e.g., by the Knife-edge and/or Uniform Theory of Diffraction 
(UTD) models, are much reduced by comparison at THz, due 
to scatterer dimensions being much larger than the wave- 
length, but nonetheless, may be considered as a form of 
multipath propagation that could contribute to the overall 
signal strength [123]. 

Moreover, the scattering behaviour of THz waves from dif- 
ferent surfaces with different roughness levels is unique at 
THz owing to the short wavelength, and a good understand- 
ing of it is required for the development of future channel 
models. EM waves that are directed on a smooth surface 
should have an equal angle of incidence and reflection (spec- 
ular reflection). However, when the surface is rather rough, 
diffuse scattering of the EM waves happens resulting in a loss 
of intensity in the specular direction. As the THz wave has a 
very short wavelength, scattering of the waves off a rough sur- 
face will cause multipath propagation and clustering effects. 
These rough surfaces can be quantified for THz frequencies, 


as having edges that vary between them less than | mm for 
frequencies > 300 GHz. Otherwise, for edges greater than 
1 mm the surface may be considered as smooth. Additionally, 
multi-layer reflections within in-homogeneous materials are 
also another possibility at THz resulting in attenuation dips 
to the received signal at selective carrier frequencies, and 
additional diffuse components [124]. 


H. Summary of Propagation Phenomena 


Table II] summarizes the approximate propagation losses at 
the center frequency of possible atmospheric windows of oper- 
ation, resulting from LoS conditions based on expected ranges 
taken from real environments for the 6G applications identi- 
fied. The propagation phenomena expected to be experienced 
by each application individually are marked by a dark circle. 
For the case of rain/snow and cloud/fog attenuation, worst case 
conditions are considered. The vast availability of the THz 
spectrum, means that certain applications can fit into specific 
windows in order to maximise the spectrum usage. For exam- 
ple, Kiosk downloads and on-board links, would suffer from 
minimal losses of 2 dB or less across the entire spectrum and 
could possibly operate in higher frequency bands. Likewise, 
intra-device links suffer from even lower losses which are 
almost insignificant. On the other hand, applications such as 
data centre, HST, V2X and WLAN communications suffer 
from moderate losses in certain atmospheric windows with the 
worst case scenario of about 57.2 dB attenuation in the W13 
atmospheric window at 30 m. Therefore these applications 
could also be utilized within a number of windows depend- 
ing on the link budget. Backhaul/Fronthaul applications, since 
they require the longest communication distances, would best 
work in the first few atmospheric windows identified. In terms 
of propagation losses resulting from scattering from objects 
found in the propagation environment, since this gives rise to 
complex scenarios that require thorough investigation before 
deployment, they are not considered in Table III. The losses 
that arise from propagation inside the human body are also 
not included for nano-network applications, since no standard 
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currently exists. With regards to attenuation due to rain and 
fog, the losses are included up to 1000 GHz for which ITU 
standards exist. Nonetheless, in the windows above 1000 GHz, 
for the applications that are expected to be utilized in outdoor 
environments, the expected ranges and corresponding losses, 
signify that these applications may not be placed within these 
bands. There is therefore a need to explore channel propa- 
gation in the windows identified, to best manage spectrum 
allocation, and maximise the capacity. 


V. THZ CHANNEL SOUNDING TECHNIQUES 


As a requirement to communication system design, pro- 
found knowledge of the propagation environment is essential. 
In order to develop future channel models, channel mea- 
surements are the first and most important step, requiring 
state of the art sounding equipment and correct methods 
to be available. THz measurements have been mainly con- 
ducted by three distinct channel sounding methods. These 
include the well known THz-TDS [125], Vector Network 
Analyser (VNA) [126] and Sliding Correlator (SC) chan- 
nel sounder [127], while THz-FDS (Frequency-Domain 
Spectroscopy) [128] is another very promising channel sound- 
ing option. In this section, the performance of these channel 
sounding methods is described and compared with measure- 
ment requirements that may vary depending on the application 
scenario under investigation. 


A. THz-TDS Measurements 


THz-TDS [125], [129], as depicted by Fig. 7(a), is a chan- 
nel sounding technique where material properties are being 
probed by ultra-short light pulses. These pulses are con- 
verted to the THz range using photo-conductive switches 
(GaAs or InGaAs/InP) irradiated by femtosecond lasers at the 
source of the Tx, and the received optical impulses detected 
at the Rx are converted to an electrical signal [130]. The 
received time domain signal is then processed by Fast Fourier 
Transform (FFT) to the frequency domain, in order to extract 
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Experimental setup for THz signal generation. (a) THz-TDS, (b) THz-FDS. 


its spectral content. Moreover, since the received signal is com- 
plex (i.e., containing both amplitude and phase information), 
it can be used to calculate the complex refractive index whose 
real part is used for measuring material thicknesses and the 
imaginary part used for calculating resonant absorption and 
extracting conductivity directly without the need for additional 
modeling [125]. THz-TDS has the advantage of large band- 
width measurements where bandwidths up to 10 THz have 
been achieved [131] allowing for a wide range of spectral 
fingerprints of different dielectric materials to be identified, 
supplemented by a very good dynamic range (~ 100 dB) 
and tuneable timing resolution that is achieved by adjusting 
the distance the two pulses will propagate; hence, enabling 
paths to be resolved within femtoseconds [125]. Additionally, 
since this is a pulsed technique, it allows for time-gating that 
can increase the Signal-to-Noise (SNR) ratio which normally 
decreases with thickness or number of layers, and allows for 
content extraction between multi-layered structures of sub- 
mm thicknesses; thus, proving to be excellent at measuring 
the scattering parameters and electrical properties of various 
material samples [132]. 

The THz-TDS poses as an excellent measurement method 
whereby characterisation for 6G applications could benefit 
from in propagation environments that are enclosed or sur- 
rounded by scatterers. The high accuracy offered by THz-TDS 
could help extract different material properties better than any 
other channel sounding technique, which can be later fed into 
channel models to characterize the propagation environment. 
It has been employed in [133] to model the reflection and scat- 
tering losses resulting from different body parts of a vehicle 
at 300 GHz for various Vehicle-to-Vehicle (V2V) commu- 
nication scenarios. THz-TDS measurements have also been 
utilized for characterizing rough material surfaces of building 
materials usually found in indoor environments, in the bands 
between 0.1-4 THz [31], [32]. Piesiewicz et al. [30] conducted 
THz-TDS measurements for future pico-cell THz indoor com- 
munication systems operating at 300 GHz. The purpose of 
this measurement campaign was to model the characteristics 
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of various building materials commonly found in indoor envi- 
ronments and the effects diffuse scattering has on the received 
power levels. 


B. THz-FDS Measurements 


THz-FDS [128] is a channel sounding technique that 
employs CW THz radiation to probe the material properties. 
The THz-FDS measurement setup is illustrated in Fig. 7(b). 
It is a technique that utilizes optical heterodyne detection in 
high-bandwidth photo-conductors that involves optical signal 
and oscillating wave mixing, where the combined product 
gives an electrical signal at the difference frequency of the 
two lasers. It offers an excellent dynamic range (~ 100 dB) 
for large bandwidth measurements (~ 6 THz). Two lasers at 
adjacent frequencies are used to illuminate the photo-mixer 
and the applied bias creates a photo-current that is modu- 
lated at the beat frequency. This allows for the antenna to 
transmit the THz wave and produce a monochromatic beam. 
Furthermore, by tuning the lasers, the THz wavelength can 
be adjusted. The main advantage of THz-FDS compared to 
THz-TDS is that is offers the user a very high frequency res- 
olution (10 MHz typ.) at a user determined range or fixed 
frequency and a comparatively lower cost [128]. Nevertheless, 
the THz-FDS technique as with THz-TDS is limited by its 
very small output power and large size and cannot be used 
for long range channel characterisation and for a wider range 
of scenarios. Furthermore, this technique cannot capture the 
dynamic behaviour of the channel such as Doppler spread, due 
to moving scatterers blocking the direct path, as it requires 
the channel to be stationary for the whole duration of the 
measurement, which may take several minutes. 


C. VNA Measurements 


The main technological advancements related to captur- 
ing the characteristics of the propagation channel are namely, 
VNA [126] and SC [134] channel sounders. The VNA, whose 
operation is more well known, works by sending a stimulus 
signal from the Tx to the Rx and processing the received sig- 
nal to obtain the channel properties. It measures the frequency 
response, which is essentially a measure of the phase and 
magnitude of the received signal by processing sequentially 
a great number of Channel Transfer Functions (CTFs) as a 
function of increasing frequency within the band of interest, 
and then aggregating them together to generate a broadband 
CTF, i.e., the S2} parameter which is a measure of the signal 
coming out of port 2 in relation to the signal entering port 
1 of the VNA and is used to characterize distance-varying 
measurements in terms of amplitude and phase, hence, the 
word vector. Onward, using Inverse Fast Fourier Transform 
(IFFT) the Channel Impulse Responses (CIRs) are generated 
as a function of time delay. A VNA can provide full synchro- 
nization between the received and transmitted signal enabling 
individual calibration and coherent averaging over long times 
for the suppression of noise because of the large number of 
narrowband measurements [135]; hence, making it ideal for 
narrowband channel characterization. To extend the frequency 
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range of the VNA to the THz range, frequency extensions are 
employed for the up-conversion of the signal [136]. 

A VNA based measurement setup is depicted in Fig. 8. 
The setup consists of a VNA connected at the Tx and Rx 
sides to frequency extender modules which in turn are con- 
nected to a waveguide and horn antenna that must be fixed 
precisely to avoid causing any signal reflections back to the 
frequency extender. To realize such a setup practically and 
enable an accurate measurement campaign and account for 
system errors in transmission and also due to source mis- 
matches, system directivity and reflection tracking [137], the 
VNA/ extender heads require lengthy calibration techniques, 
e.g., Short, Offset-short, Load and Thru (SOLT) [138]. One 
of the key limitations related to measurements using VNAs is 
that the frequency extenders have to be connected via phys- 
ical connections and require a reference signal to maintain 
phase synchronization; thus, the measurements are constrained 
in terms of Tx-Rx separation distance, while for these mea- 
surements the VNA must be placed somewhere near the Tx 
and Rx ends, therefore, acting as an unwanted reflector and 
scatterer. Furthermore, the cabled connections between Tx and 
Rx are very sensitive and may affect the calibration if they 
are moved around. Another key limitation of the VNA is the 
ability to measure the propagation characteristics over ultra- 
wide frequency ranges of a time-variant channel influenced by 
moving scatterers i.e., humans that cause blockage and scat- 
tering of the signal. In such complex cases, the applicability 
of VNAs is no longer viable since they take too long to per- 
form a full frequency sweep surpassing the channel coherence 
time [134], and therefore not allowing for dynamic channel 
variations such as Doppler shifts to be captured in order to 
analyze the effects of link blockage. Nevertheless a VNA 
can provide an excellent dynamic range (~ 120 dB) [138] 
and a very high delay resolution [124], allowing for a very 
large number of MPCs to be resolved over significantly broad 
bandwidths; Nonetheless, the bandwidth is directly dependent 
upon the selected waveguide interface [139] used for Single- 
Input Single-Output (SISO) measurements. State-of-the-art 
frequency extender modules commercially available, can oper- 
ate between 50 GHz to 1.5 THz with varying bandwidths that 
reach up to 350 GHz in the higher end [140]. 

VNAs, have been employed in a number of measure- 
ment campaigns that have been carried out for characterizing 
the propagation environment for various 6G applications due 
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TABLE IV 
PROPERTY COMPARISON OF THE MEASUREMENT CAMPAIGNS CONDUCTED FOR 6G SCENARIOS USING VNA CHANNEL SOUNDING 


Ref. Frequency Range Wavelength (4) Doppler Sampling | Delay Reso- | Path Reso- | Antenna HPBW 
Rate (< 4/2) lution lution Gain 
[141], [142] [143] | 220-340 GHz 1.071 mm 0.5355 mm 8.33 ps 2.5 mm 25 dBi 10° 
[144] 280-320 GHz 1 mm 0.5 mm 25 ps 7.5 mm 26 dBi 10° 
[137] 500-750 GHz 0.48 mm 0.24 mm 4 ps 1.2 mm 26 dBi 10° 
[145], [146], [22] | 300-312 GHz 0.98 mm 0.49 mm 83 ps 25 mm 22-24 dBi 10° 
[147] 300-320 GHz 0.97 mm 0.48 mm 50 ps 15 mm 22-24 dBi 10° 
[148] 280-320 GHz 1 mm 0.5 mm 25 ps 7.5 mm 22-23 dBi 10° 
[23] 300-312 GHz 0.98 mm 0.49 mm 83 ps 25 mm 22-23 dBi 12° 
[79] 345-355, 645-655 GHz | 0.86, 0.46 mm 0.43, 0.23 mm 0.1 ns 30 mm 25 dBi 12° 
[149] 295-305 GHz 1 mm 0.5 mm 0.1 ns 30 mm 26 dBi 10° 


to the advantages it holds of high dynamic range, high 
delay resolutions and long distance communication links. 
Measurements have been carried out for characterising the 
propagation channel for the indoor short-range kiosk down- 
loading scenario [141], [142], [143] using extender heads 
that operate between 220-340 GHz with 25 dBi horn anten- 
nas and a Half-Power Beam Width (HPBW) of 10° as well 
as for the on-board desktop scenario [144], [137] between 
280-320 GHz considering a 20 GHz bandwidth and between 
500-750 GHz using 26 dBi horn antennas with a HPBW 
equal to 10° for both cases, respectively. Several measure- 
ment campaigns have also been conducted to characterize 
data centre links using a VNA with extender heads operat- 
ing between 300-312 GHz for [145], [146], [22] and 300-320 
GHz for [147], with horn antenna gains that vary between 
22-24 dBi with a 10° HPBW in azimuth and elevation, 
for all cases. Similarly, measurements have been conducted 
for intra-device communication scenarios with an operational 
bandwidth between 280-320 GHz, but considering in the anal- 
ysis the ranges between 300-314 GHz for [148] and 300-312 
GHz for [23] with an antenna gain between 22-23 dBi across 
the frequency ranges investigated, and with a 10° and 12° 
HPBW, for each case respectively. A measurement campaign 
has been carried out for an indoor office room scenario in the 
350 and 650 GHz bands [79] with a bandwidth of 10 GHz, 25 
dBi gain horn antennas and a 12° HPBW for both frequency 
bands. An office room scenario was investigated in [149] at 
300 GHz with a 10 GHz bandwidth, 26 dBi gain antennas and 
a 10° HPBW where with the addition of lenses the gain was 
increased by 13 dBi on each side and the HPBW was reduced 
to 0.1°. The properties of the measurement campaigns con- 
ducted in relation to the frequency range of operation and 
antenna characteristics are summarized in Table IV. 


D. Sliding Correlator Measurements 


The most popular channel sounding method that can 
yield real-time results is the SC based measurement method 
that utilizes the Direct Sequence Spread Spectrum (DSSS) 
technique that was first introduced by “Cox” of the Bell 
Laboratories [150]. The SC based time-domain channel sound- 
ing method, operates by transmitting a broadband Pseudo- 
Noise (PN) sequence Srx(t) along the channel h(t) from the 
Tx to the Rx as a stimulus signal resulting to the received 
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Schematic of a Correlation based channel sounder. 


Fig. 9. 


signal b(t) and calculated by convolution of the propagat- 
ing signal with the time-varying channel. At the receiving 
side, cross-correlation of the received signal with an identi- 
cal but delayed by 7 complex conjugate of the PN sequence 
Sk,,(t — T) takes place. Given the autocorrelation function 
P(r — ¢) of the PN-sequence can be described by a Dirac- 
impulse shape which then equals to the CIR, a maximum 
length sequence (m-sequence) can therefore be used as a PN 
sequence to be transmitted as a reference signal in order 
to approximate the Dirac-impulse shaped function [60]. The 
advantage of this channel sounding method is that it can cap- 
ture all the frequency points at the same time so that it is not 
as slow as frequency scanning measurements conducted using 
a VNA. 

In Fig. 9 a simplified schematic of a SC based channel 
sounder is presented. The PN sequence is mixed up to achieve 
the targeted frequency range before modulation, known as the 
baseband. In order to estimate the complex CIRs in the base- 
band, the Tx and Rx sides need to be in synchronization 
to achieve phase recovery at the Rx side, something that is 
achieved by the clock source that is used to drive the Local 
Oscillator (LO) on the Tx and Rx sides with a high degree of 
precision. Hence, an SC can measure a dynamic channel but 
not a mobile channel due to the clock source. 

Time-domain SC operates at much higher speeds and allows 
for dynamic variations of the channel to be captured caused, 
i.e., by moving scatterers. However, strong thermal noise con- 
sequently may affect the precision in obtaining the results 
due to the lack of measurement bandwidth that is limited by 
the semiconductor technology [134]. SC channel sounding in 
the mmWave bands below 100 GHz with real-time capabil- 
ities that allow for fast electronic steering of the THz beam 
through the use of phased antenna arrays has been recently 
demonstrated [151], [152], [153]. Nevertheless, the first SC 
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channel sounder at THz has been introduced 300 GHz at 
an ultra-broad bandwidth of 8 GHz by employing frequency 
extension modules to achieve the frequency range of oper- 
ation, being able to capture 17590 CIRs per second, while 
also offering MIMO capabilities [127]. This option however, 
due to the lack of back-to-back calibration, leads to addi- 
tional unwanted MPCs, that are caused by non-linearities in the 
frequency extenders. Numerous SC measurements conducted 
by transmitting m-sequences have been recorded in literature 
using the SC channel sounder presented in [127], to sound the 
propagation channel between 302.2-308.8 GHz for complex 
environments where a large number of measurement data is 
required. Various measurement campaigns to investigate chan- 
nel properties for various HST communication scenarios have 
been carried out [154], [155], [156], [157], utilizing Tx and Rx 
antennas with a gain of 15 dBi and a 30° HPBW, while consid- 
ering varying Tx-Rx separation distances < 30 m. Similarly, 
measurements have been conducted in a real data centre envi- 
ronment with Tx-Rx separation distances < 20 m while the 
SC channel sounder has also been employed in measurements 
for various V2X scenarios [158] in a 4x4 MIMO configura- 
tion with antenna gains of 26.4 dBi and 8.5° HPBW, while 
considering Tx-Rx separation distances up to ~15 m. 

Furthermore, different sounder architectures will have to be 
developed and compared to improve the estimation accuracy 
and enhance the bandwidth. Zadoff-Chu sequences [159] may 
be used as synchronization sequences, replacing the conven- 
tional PN-sequences due to their advantageous autocorrelation 
properties. A first channel sounder has been developed at THz 
frequencies that transmits the Zadoff-Chu sequence within a 
2 GHz bandwidth over a scaleabe frequency range between 
270-320 GHz [160]. This architecture has shown to dimin- 
ish spurious components that appear above the noise floor 
caused by non-linearities experienced in [127], and has demon- 
strated a dynamic range of 70 dB, that pertains a performance 
improvement of 10 dB. SC measurements by transmitting 
Zadoff-Chu sequences have been performed using the chan- 
nel sounder presented in [160] for an indoor conference room 
scenario and an Urban street scenario [161], [162] considering 
Tx-Rx LoS distances of < 10 m up to ~ 35 m, respectively. 
Nevertheless, any results obtained using this channel sounder 
architecture cannot be considered as conclusive for accurate 
channel characterisation as they suffer from a quite low delay 
resolution of 0.5 ns that corresponds to a path difference of 
15 cm between multipaths. In addition, the acquisition time 
of 50 ms required to store 50,000 CIRs over a 1 us duration, 
corresponds to bad Doppler capabilities. A possible improve- 
ment to this system can be supported, where it could benefit 
by measuring in ‘burst, by slowing down the frequency at 
which it takes CIRs, and thus essentially reducing its sam- 
pling to handle lower Doppler shifts, to then measure for a 
longer duration. This would in turn increase its scaling factor 
which equals to 1. The capabilities of the aforementioned SC 
channel sounders are summarized in Table V. 


E. Summary 


A number of measurements have been conducted so 
far in the THz band using the above channel sounding 
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TABLE V 
PERFORMANCE COMPARISON OF THE TWO SC CHANNEL SOUNDERS 


Ref. [1] [2] 

Clock Frequency 9.22 GHz 10 MHz 
Frequency Range ~ 300-308 GHz | 299-301 GHz 
Wavelength ~ 1 mm 1 mm 
Doppler Sampling Rate | 8.8 KHz 0.025 Hz 
Delay Resolution 108.5 ps 0.5 ns 

Path Resolution 32.5 mm 15 cm 
Dynamic Range 60 dB 70 dB 
Scaling Factor 128 1 

Sequence Type PN Zadoff-Chu 
Sequence Length 4095 2000 


methodologies. Most of the measurements used to obtain 
the channels’ characteristics, have revolved around the 300 
GHz band reaching up to 750 GHz and have been con- 
ducted with the use of VNAs and SC channel sounders. On 
the other hand, spectroscopic measurements have been con- 
ducted on the entirety of the defined bands of interest between 
0.3-3 THz, mainly used to obtain characteristics of materials 
and to describe propagation at a molecular level. 

Spectroscopic measurements due to the low output power 
are suitable in characterizing short distance links but with a 
very high accuracy, and offer a huge bandwidth. They have 
been employed in characterizing various materials present in 
the propagation environment for various applications such 
as V2X and indoor communications. They have also been 
employed to obtain coefficients that are used to model attenua- 
tion by meteorological phenomena as well as for modeling the 
absorption of different tissue layers for nano-network in-body 
communications. 

A number of VNA based measurements been employed in 
obtaining channel characteristics in indoor scenarios such as 
office rooms, kiosk downloads, data centres and intra-device 
links. VNAs hold the advantage of high delay resolutions 
due to the wide bandwidths offered, allowing for reason- 
ably accurate characterisation, while the high dynamic range 
demonstrates the ability to capture multipaths over long prop- 
agation distances at scaleable frequency ranges, that would 
otherwise not be observable using the other channel sound- 
ing options. However, frequency scanning is an issue, as it is 
very time consuming. Given the wide bandwidth availability of 
commercially available frequency extension modules operat- 
ing well above 1 THz, as well as the limitations in frequency 
scaleability and bandwidth offered by SC channel sounders, 
VNAs are expected to form the backbone of future measure- 
ment campaigns that aim to obtain the propagation channels’ 
characteristics above 300 GHz. 

At THz, current SC channel sounding technologies are lim- 
ited to operating at frequencies around 300 GHz with a few 
GHz of bandwidth. These channel sounders are ideal, in that 
they can collect a huge number of measurement data in com- 
parison to VNAs, in order to model scatterer mobility and 
complex propagation environments with sufficient accuracy. 
SC channel sounders can therefore be employed to extract 
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TABLE VI 
PERFORMANCE COMPARISON OF THE FOUR DIFFERENT CHANNEL SOUNDING TECHNIQUES 


TDS/ FDS [128] VNA [126] SC [127] 
Measurement Domain Time / Frequency Frequency Time 
Bandwidth > 5 THz a few 100s of GHz 8 GHz 
Dynamic range ~ 100 dB ~ 120 dB ~ 60 dB 
Speed milliseconds to 1 min / several mins several minutes per sweep 1,000s of CIRs generated per sec 
Frequency resolution 5 GHz typ. / 10 MHz typ. 10 MHz typ. — 
Real-time Yes / No No Yes 
Range <lm a few 10s of m a few 10s of m 
Spectral selectivity No / Yes Yes No 


Intra-device links, scattering properties; 
Use Cases/Applications . E Prop 
molecular attenuation 


the channel properties for applications such as V2X and HST 
communications where the propagation environment is larger 
and much more complex than applications characterized by 
VNA measurements. 

Table VI summarizes and compares the capabilities of the 
four main channel sounding techniques described in this sec- 
tion and names a few of the applications that these methods 
would work best for in future characterisation efforts for 6G 
applications in the bands above 300 GHz. Ideally, channel 
sounders in the future would maximize these capabilities, and 
would adapt to a wider range of scenarios and potential appli- 
cations in the THz band to aid in accurate characterization and 
modeling of the propagation channel. 


VI. THZ CHANNEL MODELING METHODS 


Deriving the channel characteristics by appropriate modeling 
is the next step in communication systems design after con- 
ducting measurements. Current modeling approaches can be 
either described as deterministic [163] or stochastic, in the 
form of Geometry-Based Stochastic Models (GBSMs) [164], 
Non-Geometry-Based Stochastic Models (NGSMs) [165], 
and Correlation-Based Stochastic Models (CBSMs) [166]. 
Deterministic models, are ones that use the EM laws of wave 
propagation to determine the signal power received at a par- 
ticular location in a 3D mapped environment. Deterministic 
channel modeling approaches in the THz band involve accurate 
geometric modeling of the waves’ propagation environment and 
require exact knowledge of the boundary conditions. Accurate 
positions of the Tx and Rx antennas along with geographical 
positions and properties of different material present such as 
shape and dielectric constant are needed, to ensure the mea- 
sured and expected results are in conjunction. Even though 
they can accurately predict the characteristics of the prop- 
agating environment, deterministic models suffer from high 
computational complexity and large simulation times. On the 
other hand, stochastic models, are ones that use measurement 
data to give a statistical description of the propagation channel. 
NGBSMs are those where the scatterers are not modeled explic- 
itly and the propagation channel is considered as a number 
of delay taps, where by deriving probability density functions 
through extensive empirical data collection we can establish the 
channel model. In GBSMs probability distributions are used 
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V2X, HST communications, 
dynamic link blockage 


to describe the properties of scatterers and then fundamental 
laws of wave propagation are used to characterize the prop- 
agation process. The delay taps in NGBSMs are replaced by 
clusters in GBSMs that represent a collection of MPCs with 
similar properties, where the clusters are distinguishable from 
one another. CBSMs on the other hand, are applicable to cases 
with wide scattering, which is not evident with THz links using 
highly directive beams. Therefore, they will not need further 
consideration here. 


A. Deterministic Channel Modeling 


The most popular deterministic channel modeling approach 
that has been employed in the IEEE 802.15.3d standard, is 
known as Ray-Tracing (RT) [167]. RT in the THz spectrum 
can be employed to model various static scenarios related to 
wireless communications within, e.g., indoor office rooms, 
data centres and nano-networks. It is a method by which 
rays are launched at regular time intervals and are detected 
based on propagation in free space, where effects such as 
reflection, scattering as well as edge diffraction and cluster- 
ing are accounted for, as demonstrated in Fig. 10. Recent 
works [168], [169], have shown that RT which has also been 
employed in the optical bands, is sufficient to model and 
characterize basic propagation mechanisms in the THz bands, 
and is therefore highly applicable to simulate specific PtP 
application scenarios for operation in 6G networks at THz 
frequencies. Advanced acceleration and high efficiency 3D 
RT techniques need to be developed that can increase the 
simulator’s performance and which consider the propagation 
environments’ characteristics precisely. 

To apply ray tracing based on the environment’s geometry, 
becomes very tedious for complex scenarios and requires great 
accuracy. A strategy adopted in RT algorithms, is one which 
takes the form of a “visibility tree” [170], that holds a layered 
based structure composed of nodes and branches, where each 
node represents a scatterer, with the first node being the Tx 
and each branch represents a LoS link between the scatterers 
when that is available. When the Rx is reached, the end node 
is instead represented by a leaf. A visibility tree, can therefore 
be used to capture individual propagation paths, indicated by 
the total number of leaves created in the visibility tree, but 
can be very computationally complex. Once the visibility tree 
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Fig. 10. The RT method for deterministic channel modeling. 


takes shape, by traversing backwards from the leaf to the root 
node (Tx), and by use of geometric optic rules of propaga- 
tion including computation of the path losses, reflections using 
the Fresnel equations and accounting for diffraction effects 
by including coefficients captured using, e.g., UTD at each 
connection, the individual propagation paths can be resolved. 

An alternative approach to that of RT, is that of ray- 
launching [171], where rays are launched on a coordinate 
grid based on angular directions, following basic propaga- 
tion laws similar to RT. Simulations are therefore performed, 
where each ray path toward the grid is traced, and consid- 
ers change of direction based on a pre-set power threshold 
level within the grid, indicating whether the incoming ray has 
diverted from its path. Ray-launching can be used to sim- 
ulate practical systems with the advantage of allowing for 
multiple Rx locations to be simulated based on a single Tx, 
without laying a great deal of additional computational bur- 
den. Another deterministic approach, is that of RT map-based 
modeling such as the 5G standardised map-based models 
METIS [172], 3GPP NR channel model [173], and IMT- 
2020 channel model [174] that can generate realistic modeling 
for various application scenarios. It offers a geometrically 
simplified description of the environment, where objects are 
assigned to represent scatterers, and point sources are assigned 
to represent e.g., Tx/Rx locations, after the geometric environ- 
ment is defined, with the disadvantage of requiring knowl- 
edge of the network layout, which might not be always 
accessible. 

On the other hand, Finite-Difference Time-Domain (FDTD) 
[175] and Method of Moments (MoM), are numerical methods 


that can be used to account for back-scattering of the rays. 
They are methods that can solve Maxwell’s equations precisely 
by integral (FDTD) and differential (MoM) equation analysis. 
They can provide accurate estimates of spatial distributions 
when detailed geometric knowledge is required in terms of 
roughness of materials, shape etc. However, The amount of 
computational resources required for simulations are exten- 
sively large and the processes are time consuming, therefore 
limiting the size of the environment to be modeled. This is 
because the solution to Maxwell’s equations requires a set of 
points spaced out in the order of a wavelength so that a numer- 
ical solution to the equations can be found. This means that 
the computational complexity of the problem increases mas- 
sively as the frequency and amount of scatterers increases. 
Both the FDTD and MoM approaches mentioned are very 
attractive and may be combined with RT to provide a more 
reliable estimation method. Although, combining the two mod- 
els together may be quite challenging since they are devised 
using different tools and assumptions. Hybrid FDTD and RT 
modeling has been demonstrated in the sub-6 GHz bands [176] 
in comparison to RT only operation, showing a good improve- 
ment in terms of RMS error as well as improvements in CPU 
time. It has also been employed in the mmWave bands [177] 
by embedding diffuse surface scattering interactions obtained 
via FDTD into a ray tracer to achieve a greater accuracy in 
the results; A comparison between RT and FDTD full-wave 
simulations for intra-device scenarios at THz can be found 
in [178], where it is apparent that by just employing the 
FDTD method, simulation times would be extensively large 
even for modeling small environments, requiring a run time 
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of several months, thus, making it impractical to be applied by 
itself. 


B. Stochastic Channel Modeling 


Stochastic channel models, are ones that use empirical data 
to provide a statistical analysis of the channel by a spatio- 
temporal approach. These models are based on measurement 
campaigns and require much less complexity in terms of 
computing resource requirements compared to deterministic 
models. Although, measurement time to obtain good accu- 
racy in the results is a limiting factor. Due to the higher 
order frequencies and larger bandwidths, wideband channels 
in the THz regime require multipaths to be resolved in the 
time domain. Narrow-band fading models such as Rayleigh 
and Ricean may therefore no longer be required and a statis- 
tical model that can generate the CIRs needs to be developed. 
Therefore, in order to model the wideband channel charac- 
teristics, this would require us to obtain the spatio-temporal 
characteristics of the propagation channel. 

It can be considered that transmission is made from point 
to point and under the Wide-Sense Stationary Uncorrelated 
Scattering (WSSUS) assumption that was introduced by 
“Bello” in 1963 [95], that considers independently generated 
delay taps. A generalised form of the Bello functions can be 
seen in Fig. 11. The system function can be described as the 
Time-variant transfer function T (f, t), which is essentially a 
snapshot of the power distribution over the frequency points, 
as obtained directly from measurement. By taking the Fourier 
Transform of T (f, t), the Input delay spread function h (t, T) 
can be obtained, and by taking the spatial average over a 
local area, the Power Delay Profile (PDP) can be computed, 
which records power received over ToA of the different arriv- 
ing MPCs, where each component represents a single channel 
delay tap. The length of the PDP can be calculated as the 
number of frequency points divided by the bandwidth. Key 
parameters that can be extracted from a PDP, include time 
dispersion parameters in the form of RMS delay spread that 
represents the average received power against time delay in 
relation to the first arriving MPC, and coherence bandwidth 
that describes the frequency range over which the propaga- 
tion channel is considered to be flat. Another parameter that 
can be extracted from a PDP, is the Ricean K-factor, that can 
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Fig. 12. SV model representing exponential decay of clusters and rays over 
time. 


be defined as the ratio of the most significant MPCs’ aver- 
age power and the average received power of all remaining 
MPCs. If h (t, T) is then Fourier transformed with respect to 
time, the Delay Doppler-spread function S (7, v) is generated, 
which allows for the Doppler spectra for each of the delay 
taps to be obtained. This is useful when analysing various 
AoA/AoD and their associated powers obtained from a Power 
Angle Profile (PAP) at different time shifts, that would indicate 
which scatterers in the environment affect the wave propaga- 
tion significantly. Root Mean Square (RMS) angular spread 
is a spatial parameter that can be extracted from the PAP 
which describes the power arrival dispersed over the angles 
of incidence the antenna is being varied. Finally by Fourier 
transforming either S (T, v) or T (f, t) with respect to 7 or t, 
respectively, this then gives the output Doppler-spread func- 
tion H (f, v) that demonstrates the Doppler spectrum related 
to each frequency component. The parameters extracted from 
the above functions, are used to provide information for the 
design of transmitters and receivers that can overcome the 
effects of ISI, fading and aid in equipment calibration as well 
as in equalization processes. 

Furthermore, clustering behaviour is observed when the 
arriving MPCs or otherwise intra-cluster components have 
similar distributions in terms of ToA and AoA/AoD. Saleh 
and Valenzuela [165] proposed a widely acceptable NGBSM 
known as the Saleh-Valenzuela (SV) model that provides 
a cluster based impulse response based on reflections from 
different surfaces, where MPCs are modeled as delay taps. 
It assumes that MPCs vary uniformly in phase and have 
powers that are Rayleigh distributed and which decay expo- 
nentially based on a Poisson distributed ToA of MPCs and 
clusters. Fig. 12 shows the SV model’s impulse response that 
represents the mean square power (39(t)) over time where 
clusters are seen to decay with a rate of I with cluster arrival 
time Tı for | = 0,1,2... and MPCs are seen to decay 
with a rate of y with MPC arrival rate T4, where k is the 
MPC arrival time for k = 0,1,2... within the /*” cluster. 
Extensions of the SV model have already been developed at 
lower frequencies [179], [180], [181] to account for cluster 
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AoA in addition to its temporal behavior. Thus, PDP and 
PAP measurements may also be used as a basis from which 
the SV model can be applied at THz to fit and verify the 
resultant clustering behavior. The SV model has been applied 
successfully in modeling the clustering behavior of indoor 
environments by also associating clustering with AoA/AoD 
in azimuth and elevation through spherical coordinates, and 
was verified by empirical measurements conducted indoors 
at 300 GHz [182]. Furthermore, a modified version of the 
SV model was employed in [183] to characterize the indoor 
THz channel in the presence of antenna subarrays, validated 
based on measurements from [184] at 300 GHz, that consid- 
ers the cluster arrival time and ray arrival time within each 
cluster as a function of the number of elements of the sub- 
array. Moreover, the modified model not only considers the 
spreading loss as the conventional approach, but also caters 
for transmission behaviors by accounting for the effects of 
molecular absorption where I and y are material dependent. 
A clustering model as a modification of the SV model was 
presented in [147] for wireless data centres, where contrary 
to the traditional SV assumption, it was proposed that since 
T cannot accurately predict MPC cluster attenuations, that it 
is further separated in lı and I% as functions of time delay, 
with a selected threshold value acting as a node from which 
T can be selected based on the distribution of MPC clusters. 

The Zwick model [14], is another NGBSM that is appli- 
cable for indoor as well as MIMO channel characterization, 
that considers physical wave propagation properties and char- 
acterizes MPCs by a transfer matrix (T;(t)), that includes the 
losses, delay (7; (t)), direction of arrival (Q r ;(t)), direction 
of departure (Qp ;(t)) and scattering properties. The Zwick 
model’s non-geometric approach, can be attributed to the use 
of a spherical coordinate system, by which the locations of 
the receivers are stochastically generated whenever the direc- 
tion of motion changes as observed by the variations in AoA 
and AoD. It allows for parameter extractions for a range of 
frequencies and environments, where the evolution of MPCs 
over time is characterized using a Birth-Death (BD) process, 
a marked Poisson process. The model has been verified in the 
mm Wave bands [185], and is therefore potentially applicable 
at THz. 

GBSM models, can be categorized based on scatterer 
arrangement, with respect to the number of dimensions of 
the antenna element configuration. Based on the environment 
characteristics, GBSMs can be divided into regular and irreg- 
ular shape GBSMs, where effective scatterers are assigned to 
regular shapes in the form of one-ring, two-ring and elliptical 
distributions and irregular shapes follow a random distribu- 
tion [186]. These schemes are the most effective and reliable 
as they are highly robust against link and network failures. 
The one-ring model, may be applicable for MIMO channel 
modeling and is particularly effective for describing envi- 
ronments related to backhaul applications where the base 
station is located at a height without any obstructions and 
the access point is surrounded by scatterers, while the two- 
ring and elliptical models become relevant in indoor cases for 
6G networks where a large number of scatterers are present 
on both sides of the communications channel. Other common 
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backhaul topologies include the star, chain and tree topologies 
which are inferior to the ring topology. The star typology’s 
inferiority, can be attributed to poor frequency reuse, since all 
the links arise from the same point, making signal interference 
far more likely to occur. While the chain topology solves the 
problems related to the star topology where all Rx sites are 
connected together via a single path, it is disadvantageous due 
to the fact that when one link fails, the adjacent link will also 
fail, resulting in a complete network failure. A tree topology 
is formed as a mixture of the chain and star topologies where 
fewer links are prone to failure and require protection. 

In order to analyze the multipath behavior within individ- 
ual clusters, a second order profile for exponentially distributed 
multipaths such as the Gaussian Mixture Model (GMM) [184], 
[187], which is a cluster oriented GBSM that can be adopted 
for approximating AoA/AoD standard deviation PDFs of dif- 
ferent materials, with X} = Xə = 0, and for identifying 
similarities between the different MPCs within a cluster or 
otherwise inter-cluster components. A generalized representa- 
tion of a cluster based GBSM is depicted in Fig. 13 that shows 
a collection of rays being launched from the Tx with similar 
AoD, reflected at the surface of a scatterer and detected at the 
Rx at some AoA and ToA. As per Eq. (12), 
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these individual rays n are summed up to form a single cluster 
with a phase at the Tx (eÌ?») and mean AoD of 0, and phase at 
the Rx (ej $n) with mean AoA of @ , where the total number of 
rays is denoted as M. Then by computing the gain of each indi- 
vidual ray at the Tx (ah) and Rx (a2...) sides as a function of 
polarization ®, a single complex number is formed represent- 
ing a single delay tap of the multipath channel, where the real 
and imaginary parts are Gaussian distributed and therefore the 
single complex number formed is Rayleigh distributed. The 
Doppler effect could also be accounted for by adding the term 
elwa to the equation where w relates to the Doppler frequency. 

In [188] a GBSM was proposed for indoor AP deployment 
at 300 GHz to address the explicit relationship between ampli- 
tude and AoA for LoS and NLoS paths of indoor channels, by 
employing a zero-mean second order GMM that relates rough 
material surface scattering interactions to cluster amplitudes 
obtained via calibrated RT simulations. The model consid- 
ers AP’s located on the ceiling to guarantee the LoS path, 
and identifies the shortest NLoS path and its power for each 
subsequent receiver locations. Similarly, in [183] a zero-mean 
second order GMM is employed to model the relationship 
between AoA and amplitude in an indoor scenario at 300 GHz 
and has been found to provide a more accurate approxima- 
tion with higher angular resolutions for the propagating wave, 
especially for smaller angles in comparison to the SV model. 
Furthermore, Tekbiyik ef al. [189], presented modeling and 
analysis of short-range THz propagation channels, by utiliz- 
ing a Gamma mixture for parameter estimation in comparison 
to GMM, based on results collected from measurements in 
the 240-300 GHz band. Gamma mixture models consider a 
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TABLE VII 
STOCHASTIC CHANNEL MODELS WITH POTENTIAL APPLICABILITY AT THZ FREQUENCIES 


Model Approach | Freq. (GHz) Base Equation Ref. 
ey — Ellik)? 
Tapped delay line model GBSM 625 P(t)= 7 [95], [137] 
Saleh-Valenzuela Model NGBSM 350 B? (t) = B2 (0, 0) Zig e77 /Te-0-T)YU (t - Ti) [165], [182] 
Zwick Model NGBSM 60 T(t, f) = Frr(f).VGr(t, f)GR(t, f). 
JS CR(t. f, QR)-Trr(t, f, Qr, QR) 
.CT (t, f, Qr)dQrdQR. [14], [185] 
Zero-mean 2”4 order GMM GBSM 300 GMM (x) = oil y +Z etl j [187], [188] 
~ 2n] y2n07 ý 
Gamma Mixture Model GBSM 312 fo(x) = ZE; Pifi(x, a1, Bi), l= 1,2....L, x > 0, pı > 0 [189] 


random variable (G) that is a mixture of (L) gamma distri- 
butions considering the shape (a ;) and scale (6) parameters 
of the th component of the mixture distribution, while con- 
sidering weights (p;) that satisfy certain conditions. GMM 
has found to describe the propagation channel as well as 
the Gamma mixture model for various propagation distances; 
however, the Gamma mixture has shown a more close relation- 
ship between the actual probability distribution and predicted 
distribution of powers. The Gamma mixture model has also 
been employed successfully in characterizing the shadowing 
effects caused by multipath scattering in intra-device link sce- 
narios by modeling the PDF of the path loss fluctuations 
in [190], considering the large bandwidths and sudden changes 
in the THz communications channel. Table VII shows exist- 
ing stochastic models presented so far, detailed mathematically 
for wider applicability and comparative analysis, and shows 
the maximum frequency for which they have found applica- 
bility in characterizing THz propagation for the envisioned 
applications. 


The first stochastic channel model for kiosk applications 
was developed in [21] as part of the IEEE TEEE80.15-16- 
0168-02-003d standard as detailed in [191], based on three 
distinct geometrically derived configurations in relation to 
the Tx and Rx with varying positions, in order to enable 
fast impulse response generation to allow for an efficient 
characterization and fair comparison in relation to angu- 
lar and temporal multipath distributions for each scenario. 
Peng and Kiirner [192], presented a stochastic channel model 
for wireless data centers at 300 GHz that is included in 
the IEEE 802.15-15-0207-003d standard [191], and validated 
through delay and angular spreads by calibrated RT sim- 
ulations from [149]. The model considers LoS and NLoS 
paths, delay and angle path loss correlations as well as phase 
and frequency dispersion. Boujnah et al. [193], proposed a 
stochastic channel model including a link quality assessment 
procedure for wireless data centre environments at 300 GHz. 
The model assumes a uniform rack placement topology, where 
nodes are placed on top of the racks and THz reflectors are 
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placed on the ceiling, to reduce the effects of blockage and 
interference, and enable connections between non-neighboring 
nodes via secondary paths. Parameters considered in the chan- 
nel model include spreading, molecular losses, reflections and 
additive noise. Furthermore, the model caters for appropri- 
ate waveform selection, pulse and symbol durations that are 
tuned to reduce the system interference and optimise the link 
quality. Simulations of the proposed topology, show that it is 
possible to deploy THz communications within a data centre, 
free of significant interferers and with enhanced connectiv- 
ity. In [194] a stochastic model for indoor radio channels at 
300 GHz is presented, that was verified by calibrated RT sim- 
ulations in an office scenario. The model developed, considers 
altogether propagation in the temporal, angular and frequency 
domains in the form of ToA, AoA/AoD, frequency dependent 
phase rotations and dispersion. The model developed allows 
for stochastic simulation generation under realistic conditions, 
and caters for multiple antenna concepts. A stochastic model 
has also been developed for 270-320 GHz based on the IEEE 
802.15-15-0166-00-003d standard [191]. The model considers 
LoS and NLoS links in a predefined configuration for board- 
to-board communications and evaluates the effects walls and 
printed circuit boards have on the received signal. 


C. Combined Stochastic-Deterministic Channel Modeling 


Now what seems to be happening at THz, is that a hybrid 
modeling approach, that is, a combination of both determin- 
istic and stochastic methods is currently being preferred. To 
capture the received power at every point in space and time, 
channel sounding at THz often makes use of mechanically 
translational and rotational antennas; though, taking a very 
large sample of measurement data is often extremely time con- 
suming and modeling the environment in a stochastic manner 
with high levels of accuracy is extremely difficult. Therefore, 
in order to extend the sampling intervals within a 3D environ- 
ment to further extract the statistical behavior of the channel, 
software simulators are employed, which by themselves often 
appear rather simplistic due to lack of knowledge of the 
propagating environment and limited computational resources. 
When empirical measurements are conducted, there are often 
different components and effects that appear, which are not 
captured by RT simulators such as dense MPCs. Thus, by 
feeding the stochastic channel generator a set of empirical 
data, it can be then calibrated by replicating this data through 
extensive simulations to enable a much larger data set as an 
extension to the limited data sets collected by measurements. 

The high frequency selectivity at THz means that the prop- 
agation channel is quasi-optical by nature. This quasi-optical 
behavior, indicates that image based RT can constructively 
predict the spatio-temporal characteristics of the wave propa- 
gation and can therefore assist in quasi-deterministic channel 
modeling [195]. The Quasi-Deterministic modeling approach, 
that was adopted as part of the MiWEBA project framework 
for backhaul and access [196] based on a predefined network 
layout, allows for stochastic cluster formation consisting of the 
most significant MPCs and provides a geometric description 
of the propagation environment, considering factors such as 
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attenuation, reflection and scattering as well as blockage and 
mobility effects. An alternative to the cluster formation can 
result from use of RT to provide a mathematical center for 
the most significant MPC without the MPC actually needing 
to be present at that specific location. 

A number of stochastic THz channel generators have so 
far been developed, that are accessible to the research com- 
munity. The most widely employed stochastic simulator is 
known as cloudRT [167], which is an open source simula- 
tor that focuses on THz 6G applications. Particularly, it is a 
3D RT high performance cloud based simulator, that is cali- 
brated based on measurements conducted for various railway 
scenarios [154], [156], V2V scenarios [197], [133] as well as 
indoor scenarios [21] [123] at 300 GHz. The calibration pro- 
cess is used to form material libraries that consider constant 
EM properties along the material, based on interactions with 
various scatterers found in the propagation environment, which 
can then be employed in different environments to emulate 
different THz application scenarios. The Terasim stochastic 
simulator [198], is built as an extension to the ns-3 network 
simulator, that offers physical and link layer implementations. 
It captures the peculiarities of the THz channel for nano-scale 
and macro-scale scenarios in terms of LoS propagation, in 
relation to the effects of molecular absorption and allows the 
users to perform network selection and spectrum allocation. It 
fails however, to address mechanisms such as NLOS propaga- 
tion, misalignment and propagation based on multiple antennas 
in relation to beamsteering, beamforming, the spherical wave 
model, and does not provide an end-to-end model. The most 
recently developed channel simulator Teramimo [199], that is 
a GUI based simulator, offers 3D end-to-end modeling for 
wideband UM-MIMO propagation channels. It caters for the 
effects of LoS, NLoS, beamsteering, spherical wave propaga- 
tion, fading due to misalignment, and molecular absorption by 
complex radiative transfer theory computations following the 
mathematically exhaustive beer-lambert law up to 10 THz with 
valid approximations up to 450 GHz. Furthermore, it incorpo- 
rates capacity calculations and is verified against measurement 
based channel models found in the literature. An automatic 
planning algorithm for backhaul links is presented in [200] that 
uses data based on the IEEE 802.15.3d standard [191]. The 
planning tool employs the ring topology GBSM, and enables 
system level simulations at 300 GHz. It is seen as a suc- 
cessor to the previous proposed common backhaul topology 
known as star topology [24] due to its superior performance. 
The tool considers LoS links under the effects of molecu- 
lar absorption and the performance is analyzed in terms of 
SINR, while it evaluates the number of links that can reach the 
required SINR for wireless operation for eventual fiber optic 
replacement. This tool however does not consider NLOS links, 
which can be seen as a further enhancement. Furthermore, 
QuaDriGa [201] is a quasi-deterministic stochastic channel 
generator that contains modules which enable the implemen- 
tation of 5G standardised stochastic models such as 3GPP as 
well as map-based models etc. Recently, it has been shown 
that it can be successfully employed in the THz bands and has 
proven to be particularly successful in characterizing various 
railway scenarios such as the T2I and inside station channels 
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with a compliant framework at 300 GHz, being validated by 
measurement data and cloudRT simulations [154]. 


D. Discussion 


In 6G communication networks, the THz band will play 
a major role and will be employed for various applications 
that have been identified, and which vary by transmission 
range and propagation environment. Different modeling tech- 
niques such as deterministic and stochastic, each possess their 
own advantages and disadvantages. Deterministic models can 
provide highly accurate results but have to endure long pro- 
cessing times and a high computational burden, while they 
don’t cater for the dynamic behavior of the channel. Examples 
as such include RT and ray launching simulators employed 
to model any application environment, whether it consists of 
indoor environments (e.g., office, conference rooms or kiosk 
downloads) or outdoor environments (e.g., HST and V2X 
applications). These methods are usually calibrated by mea- 
surement data as a hybrid stochastic form of modeling in order 
to reduce the computational resource requirements, by which 
the materials the waves interacts with are characterized before 
being used in simulations to obtain the propagation character- 
istics of the propagation environment under consideration, for 
a given application scenario. On the other hand, the FDTD 
method provides more accuracy in comparison to the other 
deterministic methods and is mainly applicable in modeling 
complex structures in wideband frequency ranges. This allows 
for better characterisation of the scatterers present in the envi- 
ronment and highly applies for indoor THz applications where 
the amount of scatterers is large and where interaction with 
them cannot be avoided, as not only the received signal is 
subject to additional propagation paths, but is expected to be 
utilized to contribute in some cases (e.g., data centre communi- 
cations) by D-NLOoS reflections. Stochastic models can provide 
a statistical description of the environment, though uncertainty 
in the developed model is obvious, due to lack of knowledge 
of the propagation environment and due to the limited chan- 
nel sounding opportunities that ensue channel sparsity as a 
result of the high directionality of the THz wave. Most of the 
stochastic models applied follow the WSSUS assumption, and 
make use of the Bello functions to analyze the channel and 
obtain its characteristics. NGBSMs such as the SV model, 
have been rigorously tested for high frequency applicability 
at THz to model the clustering behavior, while the Zwick 
models’ spherical coordinate approach allows for cluster evo- 
lution to be captured. GBSM approaches such as the GMM, 
are applicable in indoor environments to model multipaths 
within individual clusters, while recently it has been shown 
that Gamma mixtures outperform GMM. The star network 
topology GBSM seems the most robust in modeling Backhaul 
applications in comparison to chain and tree topologies and 
can be employed in simulators to generate the stochastic 
configuration and evaluate the SINR of the channel. 


VII. THZ CHANNEL CHARACTERIZATION FOR 
6G APPLICATIONS 


In this section, the channel characteristics for each of the 
application scenarios identified in the THz band are presented 
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in relation to the propagating environment under consideration, 
and the key findings are discussed. The channel character- 
istics presented, include fading characteristics such as path 
loss, and dispersion characteristics that affect system design 
such as delay spread, angular spread, coherence bandwidth 
and K-factor. 

1) Short-Range Channels: Song et al., presented a VNA 
based channel measurement campaign for the kiosk download- 
ing scenario at 240-320 GHz [143]. Initially, free space loss 
exponents were measured within the range of 1.8-2.3 for the 
entire frequency band. Due to the cavity effect between the 
transmitter and the housing package employed to emulate the 
kiosk, three late multipaths arrived that were almost eliminated 
when the housing of the receiving antenna was tilted. The BER 
was then measured against separation distance between the Tx 
and Rx, where 20 Gbps error free transmission was achieved 
over | m. He et al. [21], [141], obtained channel characteristics 
in the 220-340 GHz band for 3 different Kiosk download- 
ing scenarios, (i) LoS propagation, (ii) propagation where the 
transmitter side window is tilted and (iii) when the receiver 
and receiver window are tilted at the same angle. For scenario 
(i), three secondary reflections arrived, that contributed to ~ 
30% of the overall signal strength, along with minor fourth 
order and sixth order reflection contributions. For scenario (ii), 
many of the reflected paths where suppressed, with second 
order reflections contributing to about 8.82% of the overall 
signal strength along with negligible fourth order reflections, 
while for scenario (ii), multipath propagation was essentially 
eliminated. In terms of RMS delay spread modeling, average 
values were seen to decrease, with values of 0.39, 0.22 and 
0.12 ns, while average Rician K-factor values increased, with 
values of 14.91, 22.13 and 26.07 dB, for scenarios (i)-(iii), 
respectively. 

Kim and Zaji¢ [144] presented a statistical characteriza- 
tion for propagation channels on-board a desktop in the 
300-320 GHz band for LoS and NLoS communication scenar- 
ios. For the LoS scenario, the path loss was found to increase 
significantly with frequency without following the theoretical 
Friis prediction due to losses in the transceiver. For the NLoS 
scenario in a specular configuration, the path loss from reflec- 
tion off an FR-4 substrate and plastic material showed close 
match to theory, while reflection off a metal, strangely, showed 
an additional loss of 5 dB. For the LoS scenario, the delay 
spread followed a normal distribution over all measurement 
distances with values ranging between 4.35 to 484.94 ps, and 
coherence bandwidth values between 36.6 to 0.328 GHz. The 
results showed an increase in delay spread with an increase in 
path loss. For the NLoS scenario, delay spreads of 210, 187 
and 227 ps, as well as coherence bandwidths of 0.757, 0.85 
and 0.701 GHz were measured. In terms of multipath prop- 
agation, any components were seen to diminish by including 
absorbers. Short-range characterization was further extended 
to a carrier frequency of 625 GHz with 250 GHz bandwidth to 
investigate the channel properties in the ultra-wideband [137] 
for the on-board desktop scenario. By analyzing the PDP, three 
delay taps were seen to arrive, lagging the main component 
by ~ 1.34, ~ 1.74 and ~ 2 ns with decreased powers of about 
30 dB, which were calculated to be losses by reflections off 
the desk, the frequency extenders as well as from the VNA. 
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From the above it can be concluded, that for the Kiosk 
download scenario, the presence of a front cover window for 
LoS cases, would induce multiple order reflections with paths 
arriving after the main component. Wave propagation would 
benefit by tilt of the transmit antenna as well as tilting of the front 
cover window to suppress these reflections. Furthermore, due to 
constraints on the experimental equipment that require housing 
to emulate the kiosk system, measurements at different angles in 
azimuth and elevation have not been realized and require further 
consideration. on the other hand, communications on board a 
desktop, between devices such as laptops PCs and routers, do 
not suffer from packaging constraints but the height of the 
antennas and surrounding equipment would cause secondary 
components reflecting off the desk and surrounding scatterers, 
and must therefore be carefully positioned in close proximity 
to minimise the losses resulting from multipath scattering. 

2) Data Center Channels: Cheng et al. [145] presented 
a THz MIMO channel sounding measurement campaign in 
a data centre like environment at 300-312 GHz for the 
LoS and Obstructed-Line-of-Sight (OLoS) propagation scenar- 
ios, where cables were considered as the obstructions, with 
measurements conducted for separation distances between 
40-100 cm. Path loss exponents were measured around 1.9 
for both scenarios, shadowing gain was measured as 0.06 
and 2.6 dB, for LoS and OLoS, respectively, while delay 
spreads were seen to be log-normally distributed. Moreover, 
Cheng et al. [146], presented an amplitude fading analysis for 
a 4x4 system using a virtual uniform linear array configu- 
ration for the 300-320 GHz band for an OLoS scenario with 
stationary cables, as well as Doppler modeling for a scenario 
with vibrating cables due to the effect of cooling airflow tur- 
bulence. For the first scenario, at a Tx separation of 15 cm, 
the calculated mean path loss was measured at —87.75 dB 
and the logarithmic standard deviation of the shadowing gain 
at 6.03 dB. The RMS delay spread was calculated to have a 
mean value of —94.13 dB and standard deviation of 2.31 dB. 
For the Doppler analysis related to the second scenario, the 
max Doppler frequency observed was 56.1 Hz, corresponding 
to a scatterer velocity of 5.61 cm/s. 

In addition, Cheng et al. [147] extended their analysis in 
the 300-312 GHz band in a real data centre environment, 
with servers and racks at larger distances in comparison to 
previous campaigns for the scenarios of (1) server to server 
LoS link, (2) OLoS link between mesh door and cables, 
(3) NLoS link by reflection, (4) Obstructed-Non-Line-of-Sight 
(O-NLoS) obstructed by mesh door and server rack, and (5) 
LoS on top of a cooling grill. Measured average path losses 
and delay spreads at 175 cm were 57.7, 63.4, 78.4, 84.4 dB 
and 0.295, 0.113, 0.372, 0.299 ns for (1) door open, (2) door 
closed, (3) door open and (4) door closed, respectively. Some 
multipaths could be seen for cases (1)-(4), 20-40 dB weaker 
than the main peak. The impact of cooling on air flow was 
found to be negligible. Furthermore, a clustering analysis was 
performed, where the average Poisson distributed arrival rate 
(à) was found to be 1.64, with a mean number of clusters 
(Neluster) being 2.64. 

Most recently, Cheng and Zaji¢ [22] presented a 
characterization of propagation phenomena for the cases of 
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(i) Rack-to-Rack (R2R), and (ii) Blade-to-Blade (B2B) links 
in a wireless data centre at 300-312 GHz. The different 
scenarios considered include LoS, OLoS, Reflected-Non-Line- 
of-Sight (R-NLoS), Obstructed-Reflected-Non-Line-of-Sight 
(OR-NLoS) propagation for each of the two cases. For R2R 
LoS and R-NLOoS, path loss values for 50-210 cm distances, 
ranged between 50-64 dB. PDPs measured showed the pres- 
ence of multiple weaker reflections. For B2B R-NLoS and 
OR-NLOS cases, the path loss was seen to fluctuate as distance 
increased, while misaligned reflectors caused an additional 
path loss of 4 dB above Friis’ values. Furthermore, the B2B 
OR-NLOS link showed that cables as obstruction can increase 
path loss by around 5-10 dB as opposed to the R-NLOS case. 
R2R OLoS and R2R OR-NLoS Links showed an additional 
5-10 dB path loss in the presence of obstructing cables, while 
OR-NLOS was seen to outperform the OLOS link with 2.5 dB 
weaker shadowing gain and weaker multipath. Measured PDPs 
for B2B R-NLoS and B2B OR-NLoS at 24, 32, 42 cm 
distances showed a cluster forming at twice the main com- 
ponents delay, which for the OR-NLOoS case, diminished in 
the presence of cables. 

Despite the previously mentioned efforts, Eckhardt et al. 
[142] conducted their own measurement campaign at 300 
GHz in a real data centre. Reflections off the racks showed 
two components arriving after the main peak, being just ~ 5 
dB weaker in strength. Measurements were also conducted 
between curtains that connect the server racks to the ceil- 
ing, where the signal underwent an additional attenuation of 7 
dB. Furthermore, specular reflections were measured through 
consecutive rows of curtains were the power reflected by the 
second and the third rows decreased by 11.4 and 17.7 dB, in 
comparison to the specular reflection resulting from the first 
row, with a distance between the curtains of 4.114 m and 1.890 
m, respectively. 

From the above findings, it can be said that propagation 
within data centres can be considered as a complex scenario 
due to many peculiarities found in the propagation envi- 
ronment. The presence of cables and mesh doors acting as 
obstructions in the propagation environment cause additional 
losses. D-NLoS links can heavily contribute to the received 
power, inducing minor losses compared to the LoS links by 
a few dB, indicating that they can be employed to efficiently 
contribute to the received signal without the need for expen- 
sive technologies such as RIS. Cabled obstructions have been 
observed to reduce multipath propagation as a trade off for 
weaker received power from the main component. Finally, tur- 
bulence due to airflow causing cable vibrations, may induce 
Doppler shifts that may negatively affect post processing of 
the received signal. 

3) Intra-Device Channels: Kim and Zajié [148], presented 
channel characteristics of propagation on a computer moth- 
erboard at 300-314 GHz. Five different scenarios were con- 
sidered, (i) LoS with a large ground plane, (ii) LoS with 
Tx/Rx height variation, (iii) R-NLoS via vertically placed 
components, (iv) OLOS, i.e., between two vertically inserted 
components, and (v) NLoS via a heat sink and rotating fan. 
Path loss for (i) was found ~ 4 dB above Friis’ prediction, 
while for (ii) at 0.8 cm and 2.1 cm Tx/Rx variation, path loss 
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was similar to the theoretical value. All path losses across 
the entire frequency range were seen to fluctuate with up to 
5 dB variations. RMS delay spread values and coherence band- 
widths were recorded as 118, 137, 161 ps and 1.35, 1.16 and 
0.99 GHz for Tx/Rx height variations of 0, 0.8 and 2.1 cm 
respectively. For scenario (iii), a RAM card was considered as 
the reflector, where the highest path losses of ~ 30 dB above 
the Friis’ theoretical value were found for angle deviations of 
48° and 36° and the lowest path losses of ~ 7 dB above theory 
for angles of 40°, 42°, 44°. The PDPs for various Tx/Rx sep- 
arations showed a relationship between multipaths arriving at 
longer delays with greater height variations due to the longer 
path lengths. For scenario (iv), path loss was seen to increase 
with an increase in Tx/Rx height separation as a result of 
longer path lengths. Finally, for scenario (v) and for different 
Tx-Rx separation distances, the path loss observed was less 
than Friis’ prediction due to constructive interference between 
the LoS path and multiple bounces, besides the halfway offset 
point, where destructive interference caused the path loss to 
be above Friis’ prediction. 

Fu et al. [23], presented channel characteristics for measure- 
ments conducted within a desktop size metal enclosure (30.5 x 
30.5x 10 cm) at 300 GHz. Five distinct propagation scenar- 
ios were considered, (i) LoS, (ii) R-NLoS with and without 
memory card as the reflector, (iii) OLoS with parallel plate 
structure as the obstructor as well as (iv) NLoS with heat sink 
as the obstructor, and finally, (v) LoS over an FPGA board. For 
(i), with Tx/Rx height variations < 1.8 cm, the path loss was 
less than Friis’ theoretical value, while for variations >1.8 cm, 
they were in line with Friis’ value, with fluctuations of 3-4 
dB across the entire frequency range. The PDP for an empty 
enclosure showed several clusters of multipaths arriving within 
2.05 ns after the main component due to the signal bouncing 
around the enclosure. For scenario (ii), the path loss decreased 
by about 15 dB with the presence of the card, while the more 
rough side of the card showed 5 dB greater path loss than 
the flat side. From the PDPs, the rough side showed multi- 
paths arriving 1.17 and 2.05 ns later and for the smooth side, 
multipaths arrived 2.05 and 4.09 ns later. Delay spread values 
of 0.1251 and 0.0618 ns and coherence bandwidths of 1.273 
and 2.576 GHz were recorded, for the flat and smooth side, 
respectively. For scenario (iii), path losses for 4 cm plate sep- 
aration were less than Friis’ value, while for > 4 cm, they 
were aligned with theory. Coherence bandwidth with increas- 
ing distance separations tended to zero, while delay spreads 
increased from 0.07 to 0.4 for 2 to 7 cm distance separation 
between the plates. Five late arrival peaks were observed from 
the PDP for all plate separation distances. For scenario (iv), a 
destructive effect due to propagation within the fins of the heat 
sink and constructive effect when reflected off the fins, and the 
path loss was lower than theory due to a waveguiding effect 
in the first case and due to resonating modes in the second 
case. Finally for scenario (v), the path loss was aligned with 
Friis’ theoretical value with four delay taps occurring after the 
main peak. 

Very recently, Fu et al. [190] aimed at addressing out- 
standing issues that have not been considered in prior 
works, conducting measurements between 300-312 GHz, 
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considering the effects of circuit boards and motherboards 
placed within the metal cavity in comparison to the free- 
space case without the metal cavity. The scenarios considered, 
include (i) LoS over the motherboard, (ii) R-NLoS with 
memory card as the reflector, (iii) OLoS with two parallel 
memory cards as the obstructors, (iv) NLoS with a heat sink 
as the reflector and finally, (v) a practical CPU-PCI link. For 
scenario (i), measured path loss between O and 2.4 cm fol- 
lowed theory, except for heights of 0, 0.6 and 1.2 cm where 
additional attenuations of 2-3 dB were caused in free space as 
well as due to reflections from the motherboard and resonances 
in the presence of the cavity. For the free-space case in (i), a 
single large multipath cluster was seen, while in the presence 
of the cavity, multiple reflections caused additional MPCs. For 
scenario (ii), considering both sides of the reflecting card, the 
path losses with and without the cavity were found greater 
than Friis’ prediction, with the rough side showing the largest 
path losses. From the PDPs, two late clusters were observed 
in free-space due to reflections from back panels, while in 
the presence of the cavity, four late arriving clusters were 
observed, due to the side wall and bouncing off the cavity. For 
scenario (iii), by varying the distance between the obstructing 
plates, the path loss followed the Friis prediction, except at 2 
cm separation for both cases with and without cavity, due to 
destructive interference causing additional losses of 10-15 dB 
above theory, with heavy fluctuations in free-space, and a 
linear increase with additional losses of 3-8 dB with the cav- 
ity. The PDP in free-space showed multipaths resulting from 
bounces between the back panels and the PDP with the cavity 
showed late arriving clusters due to bounces resulting from 
the metal cavity. For scenario (iv), all measured path losses 
were above Friis’ prediction with heavy fluctuations, espe- 
cially for the case with the metal cavity. The PDP in free-space 
showed three late arriving clusters due to an extra propagation 
distance between the heat sink and back panel as well as dis- 
tance between back panels and antennas, while the PDP with 
the cavity, again showed multiple clusters arriving late due to 
the signal bouncing around the cavity. Finally, for scenario 
(v), strong fluctuations were seen on the path loss which was 
above Friis’ prediction for both cases, with additional losses 
between 15-25 dB. PDPs showed similar behavior as before, 
with two clusters arriving in free-space due to reflections from 
the back panels and multiple late clusters arriving for the case 
with the cavity, resulting from reflections from the sidewalls. 

From the above findings, it is apparent that for the case 
of propagation within a metal enclosure, additional multipaths 
would be expected resulting from reflections in comparison 
to free space, while cavity resonances may cause additional 
losses and fluctuations. NLoS links were seen to have the high- 
est path losses from the scenarios considered, while proper 
placement in terms of height of the antennas are required to 
minimise losses in LoS links, where even lower values than 
Friis’ predictions can be achieved. Propagation between par- 
allel plates in the OLoS case or in any other similar case, 
could benefit from the small wavelength as such to form a 
waveguiding effect that would ultimately increase the received 
power. Furthermore, components such as RAM cards showed 
that additional losses and multipaths can be expected and 
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require further characterisation in terms of roughness for the 
R-NLoS case. Moreover, devices such as heat sinks, as well 
as back panels contribute to additional late multipaths arriving 
for R-NLoS propagation. By proper placement of the anten- 
nas, losses can be minimized in order to eventually achieve 
extremely large data rates, given that atmospheric losses at 
these distances are negligible. Generally, due to the confined 
space and small wavelength of the THz wave, clusters of multi- 
paths were formed, and broadened peaks were observed due to 
diffuse reflections occurring from the roughness of the compo- 
nents. Lastly, it can be said, that the aforementioned results are 
adequate in characterizing intra-device channels at 300 GHz 
and can be used as a basis for characterizing intra-device 
channels in higher frequency bands. 

4) Backhaul/Fronthaul Channels: Based on current 
demonstrations, measurement distances up to 120 m can be 
realized using the RFoF technique, and thus realistically, 
the envisioned application of backhaul/fronthaul must be 
within that range. At 300 GHz, the signal traveling at LoS, 
will experience atmospheric attenuation losses of ~ 4-5 
dB at 1 km, which means that at a distance of 120 m, 
the attenuation formed will be ~ 0.5 dB. The main signal 
attenuation for this case, is the FSPL which is ~ 15 dB. 
Depending on the height of the Tx and Rx, the propagating 
wave may experience the effects of rain that can induce an 
additional attenuation of ~ 2.5 dB for extreme rain conditions 
at 300 GHz, while the effects of fog, are small enough (< 
0.1 dB) that can be neglected. 

Since measurement range capabilities are limited based on 
the current technological advancements, automatic planning 
algorithms based on network topologies are employed, that 
consider realistic propagation environments with obstructions, 
to help identify possible wireless connections in order to eval- 
uate the SINR characteristics of the propagation channel and 
achievable data rates. For example, in [200], a realistic com- 
plex scenario of a dense cellular network in Hanover, Germany 
was considered for 300 GHz, with an average link distance of 
144 m, an antenna gain of 50 dBi and varying meteorolog- 
ical conditions. The tool considers the maximum achievable 
distance of a wireless link, and assumes an LoS condition 
with a safety margin in terms of angular displacement for the 
case where the antennas are not perfectly aligned. Comparing 
the two network topologies, namely, ring and star, the former 
results in twice the number of established connections in com- 
parison to the later one. The SINR values of all the links in the 
ring topology are seen to mostly exceed 22 dB with less than 
1% failures, while 90% of the links in the star topology do 
not exceed 22 dB. These results were generated considering a 
5 GHz bandwidth and 50 GHz bandwidth for the worst case 
conditions of signal attenuation. 

Automatic planning simulation tools are extremely impor- 
tant when designing communication systems to operate for 
backhaul/fronthaul applications. It has been observed that even 
with extreme weather conditions, link stability is maintained, 
with the ring topology posing as the best choice due to the 
improved SINR of the channel, and therefore the feasibility of 
such systems is proven to be a solution for fibre-optic cable 
replacement. Furthermore, outdoor measurements at 300 GHz 
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and above are missing from literature and must be further 
explored to verify the applicability of these tools, and allow 
for a much larger propagation characterization, that would 
improve communications system design for backhaul/fronthaul 
applications. 

5) WLAN Channels: Priebe et al. [149] investigated the 
propagation channels’ characteristics at 300 GHz within a 
small indoor office room of size 2.43 x 4.15 m for vari- 
ous AoA and AoD for both the Tx and Rx, to emulate fixed 
PtP communications between a mobile device and an access 
point. For the direct LoS path, a deviation of about 3.9 dB was 
observed between measurement and theory for a 1.67 LoS link, 
while similar losses can be assumed due to misalignment for 
the D-NLoS case. For various combinations of AoA and AoD, 
corresponding to reflections off a wall, wardrobe and door, 
the PDPs demonstrated that the NLoS component provides 
insignificant contributions to the received power. The only 
multipaths that occured were for the LoS case when varied is 
small angular steps up to 5° to cater for misalignment effects, 
which nonetheless, were seen to be very weak. Additional 
losses in addition to Friis’ predictions were recorded between 
13.9 and 30 dB, as well as between 29.7 and 52.7 dB for the 
primary and secondary paths, respectively. RMS delay spread 
values were calculated, and ranged between 72.1 to 127.5 ps 
and 112.1 to 117.2 ps for the primary and secondary paths, 
respectively. For the LoS case, the RMS delay spread was 
recorded as 16.5 ps. 

Priebe and Kurner [194] provided a statistical analysis for 
an indoor propagation channel at 275-325 GHz for a realistic 
office room of size 3.59 x 45.2 x 2.82 m that consisted of 
wardrobes and tables. The scenario was intended to emulate 
the case of communications between an access point below the 
ceiling in two different locations, one in the centre and one in 
a top corner, and a user device at different locations. Double 
directional data were recorded by varying the Tx and Rx in 
azimuth and elevation. In terms of LoS propagation for first 
and second order reflections, the Gaussian distribution was 
found to fit well the relative amplitudes and amplitude dif- 
ferences, while temporal decay was approximated well with 
a Rayleigh distribution and amounted around 1 dB/ns. Ray 
deviation from a straight line was recorded between —10 and 
10 dB. First order and second order reflections were found 
to be attenuated by more than 10 dB and 20 dB on aver- 
age in comparison to the LoS path. Amplitude decay was 
recorded between 0.25 to 0.5 dB per ns while amplitude devi- 
ations ranged between —5 to 5 dB and —7.5 to 7.5 dB for the 
two different access point locations. AoA angular spread for 
azimuth and elevation, showed deviations of 10° and higher 
up to 40° depending on the Tx horizontal location and up to 
a few degrees, respectively. For the access point scenario that 
was close to the wall, higher delay spread deviations occurred. 

Zhao et al. [79], extended THz indoor channel characteriza- 
tion to higher frequencies by providing results for a conference 
room scenario of size 9.20 x 5.98 m in the 650 GHz band, 
and in comparison to the 350 GHz band. The Tx was placed 
in one corner for two Rx locations, one in the opposite cor- 
ner of the Tx and the other one on the diagonal of the Tx. 
Ten resolvable paths were found, including LoS and up to 3rd 
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order reflections were recorded for specific traceable AoA and 
AoD combinations, by considering a full 360° azimuth rota- 
tion. Propagation losses for the 650 GHz band in comparison 
to the 350 GHz band, for path lengths between 6.81 and 15.35 
m were recorded to be greater by 0.6 to 6.3 with no linear 
relation between path length and loss. Extrapolated channel 
capacity in Gbit/s was measured as being less than | for all 
resolvable paths except the LoS path for the 350 GHz band, 
and less than 0.01 for all resolvable paths except LoS and one 
first order reflection path for the 650 GHz band. PDPs were 
also recorded for the same resolvable paths as for the PAPs, 
where no significant multipaths were recorded for all resolv- 
able paths in both the 350 and 650 GHz bands, and therefore, 
the channel can be considered to be flat. RMS delay spread 
were measured as < 0.08 and < 0.06 and are seen to fol- 
low Gaussian and Rayleigh distributions for the 350 GHz and 
650 GHz bands, respectively. Coherence bandwidths were then 
calculated as < 5 GHz, except for one of the two LoS paths 
for the 350 GHz band, and equal to 3.3 GHz except for one 
LoS path and two first order reflection paths with 10 GHz 
coherence bandwidth for the 650 GHz band. 

NLoS propagation is regarded as a method by which the 
THz wave can benefit from besides contributions of the direct 
path. Nonetheless, it is not clear yet as to whether NLoS prop- 
agation can actually contribute to the received power in such 
applications, even for a small in size propagating environment. 
This comes from the fact that materials in the propagating 
environment differ in characteristics, and need to be mod- 
eled accordingly, to determine their effect in each scenario 
and environment. Materials with roughness in the order of the 
wavelength and above, will diffusely scatter the power to a 
greater extent, and therefore become more prominent in the 
higher bands, while materials with higher reflective properties 
will have better contributions to the received power by NLoS 
propagation. Contribution by NLoS propagation also depends 
upon selecting the best Tx to Rx configurations, while main- 
taining at the same time a good LoS alignment where possible. 
Furthermore, the signal should only be received within a small 
range of angles, to actually contribute when the LoS path is 
blocked. Nonetheless, D-NLoS can be considered as the most 
beneficial mechanism besides LoS that can be utilized when a 
direct link is absent, and hence a greater need to characterize 
the materials the wave interacts with. Lastly, misalignment will 
cause large losses in the power received and must be strictly 
avoided. 

6) HST Channels: Guan et al. [154], characterized the T2I 
inside-station channel for the first time at 300 GHz. Two dis- 
tinct scenarios were considered, (i) Tx close to a lamppost 
that was near to a parallel train to emulate a relay point with 
the Rx near to the wagon with the door open, and (ii), the Tx 
close to the lamppost and the Rx was placed between the two 
wagons. Based on the PDP measurements, the Rician K-factor 
and RMS delay spread were measured as 3.52, 3.60 dB and 
8.92 and 9.23 ns for (i) and (ii), respectively. These values 
demonstrate a good LoS condition and greater multipath spar- 
sity compared to the lower mmWave bands. Interestingly, for 
the first measurement, a weak reflection occured from the inner 
wagon wall through the open door arriving with a long delay 
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of more than 60 ns. Furthermore, based on extensive simula- 
tions where the Tx was placed on the platform and catenary 
mast, the path loss showed a 2-Ray model pattern due to the 
constant fluctuations when the Tx is deployed at the platform 
in (i), while the fading is much more random when the Tx is 
deployed at the catenary mast between two tracks as in (ii). A 
low K-factor with median value below 0 dB and small delay 
spread was measured with LoS, indicating that much stronger 
multipaths occur when the Tx is employed on the catenary 
mast due to reflection attenuations caused by the metallic body. 
Furthermore, the azimuth angular spreads of arrival were much 
greater than those of departure, meaning that most of the scat- 
terers are surrounding the receiver and most likely come from 
the train body. Moreover, the angular spread of departure is 
small in the catenary mast when no other train is present in 
the station. When another train was present, the angular spread 
jumped from 1° to 40°, causing additional multipaths. 

Guan et al. [156], [155], investigated the intra-wagon use 
case at 300 GHz. The scenarios considered include LoS, 
NLoS, Light-NLoS without LoS but with the first order reflec- 
tion, and Deep-NLoS for the NLoS region without the first 
order reflection. A coverage analysis was conducted based 
on extensive calibrated simulations, which showed that three 
access points will be needed to optimally enhance the coverage 
in one wagon. This comes from the fact that deploying two 
Tx’s on the side and one in the centre of the wagon in accor- 
dance with the LoS and Light-NLoS contributions, results in 
an area coverage of 78.1%, in comparison to deploying the Tx 
only in the centre or only on the side of the wagon with an area 
coverage of 49.5% and 36.9%, respectively. Ricean K-factor 
measured for the Deep-NLoS region was 8-13 dB lower than 
in the Light-NLoS region, while in the Light-NLoS region 
and LoS regions, it differed by only 2-6 dB, with the Ricean 
K-factor for the LoS region measured between 4.27-6.9 dB. 
These observations indicate that the first order reflection is 
the dominant propagation mechanism when the LoS is absent. 
Very small RMS delay spread values were measured between 
1.15-2.57 ns, that are largely influenced by the propagation 
environment, deployment strategy and frequency of operation, 
as similar sort of RMS delay spreads were measured in the 
sub-6 GHz bands [202]. In-terms of angular spreads, the Deep- 
NLOoS scenario showed much larger standard deviations than 
the LoS and Light-NLoS scenarios, indicating that propaga- 
tion heavily relies on the nature of the weaker scattered paths 
for the Deep-NLoS case. Furthermore, when multiple strong 
reflections occur, this causes the mean angular spreads for 
the Light-NLoS region to be much greater than the LoS and 
Deep-NLoS regions, meaning that THz links can be further 
enhanced when the LoS is not present. 

Guan et al. [157], presented channel characteristics and 
a capacity analysis for all of the identified railway scenar- 
ios, namely, T2I, inside-station, T2T, I2I and intra-wagon. It 
was determined that 99% of the channel capacity is much 
lower than the average channel capacity when the Tx and Rx 
are aligned, and this can be attributed to multipath fading. 
This effect requires highly directional antennas to be miti- 
gated. Furthermore, The mean Doppler shift and RMS Doppler 
spread values in the T2I and T2T channels were found to be 
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73.2 KHz and 8.25 KHz, respectively. These values indicate 
that this is beyond the capability of Doppler shift calibration 
requirements in the current 5G equipment and therefore must 
be considered carefully in the design of THz communications 
systems for railways. 

The results show that different deployments and setups 
result into different channel characteristics, and therefore 
extensive measurement campaigns are required to identify the 
best deployment strategy for each scenario. An RIS could 
be effectively employed in a D-NLOoS scenario to enhance 
the first order reflection when the LoS path is not available. 
Nonetheless, unique reflections such as from the inner wagon 
wall when the door is open during the time the train is sta- 
tionary in the T2I scenario must be considered before system 
deployment as they may effectively cause ISI. The extensive 
simulation results for all five smart railway scenarios con- 
sidered can be fed into channel generators that would bring 
the channel characteristics or channel models into the stan- 
dard channel model families like in 3GPP, which will be used 
for further evaluation of the system design and for future 
performance analysis of the communications system. 


7) V2X Channels: Petrov et al. [133], conducted an 
interference study that considered the impact of direct trans- 
mission and multipath reflections from vehicles on the side 
lane for urban and highway V2V communications scenarios 
at 300 GHz. The scenarios considered include LoS, and angu- 
lar surface scattering from reflections off (i) Aluminium for 
the cars’ main body, (ii) tempered glass for the windows, and 
(iii) plastic for bumpers and fenders, which are parts typi- 
cally found on a car. The measurements were conducted for 
angles between 35° and 75° that would in practice have the 
greatest contributions. For (i), ~ 0 dB attenuation at 35° was 
measured, that was increased to 3 dB at 75°. Equivalently, 
attenuation values for (ii) were recorded as 2.5 and 13.2 dB, 
while for (iii), the corresponding attenuations were found as 
7.3 and 14.6 dB. Furthermore, to assess dense traffic condi- 
tions, neighboring vehicles separated by 1 m from each other 
were placed at a 0.2 and 0.3 m offset to the Tx and Rx com- 
municating vehicles, that were placed on a clear LoS path and 
bumper height, and were varied from 10 to 25 m. Up to 2nd 
order reflected and scattered rays were observed, where the 
difference between them and the LoS component in terms of 
power levels were found to be ~ 5 dB and ~ 25 dB at 10 
m that increased over distance to ~ 12 dB and ~ 35 dB at 
25 m. Furthermore, it was found that the interference largely 
depends on the antennas’ directivity angle, where for angles 
greater than 40° multipath reflections from the side lanes con- 
tributed significantly to the interference, while for angles less 
than 20° the interference caused was minimal. 

Eckhardt et al. [158] conducted a detailed propagation anal- 
ysis based on measurement data obtained at 300 GHz for V2V 
communications in an indoor stationary environment. For clear 
LoS propagation, measured path losses showed an additional 
expected attenuation error < 3 dB for different height vari- 
ations, and ~ 1 dB due to alignment for measurements at 
the same height compared to Friis’ prediction. For LoS vehi- 
cle blockage, attenuations of 40 and 50 dB were recorded for 
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1 and 1.2 m height, while at the windshield level, attenua- 
tions dropped by 20 dB at 1.25 m and then increased to 45 
dB at 1.4 m. For under-vehicle propagation, at 16 m Tx-Rx 
difference, for different heights < 0.9 m, losses varied from 
64 dB to 0.5 dB at 0.4 m. For single-lane front and rear reflec- 
tions, measurements were taken at a Tx/Rx height of 58 cm 
for distances between | and 15 m. The impulse responses 
for all measurement distances showed a broadened peak that 
decreased to noise level at about 10 ns. The difference between 
the path loss and FSPL was constant, recorded at about 25 dB 
for all distances due to reflection losses. RMS delay spreads 
for the rear reflection averaged at 1.22 ns, ranging between 
0.73 to 2.44 ns. For a multi-lane scenario, the reflection loss 
averaged at 5 dB with a minimum of 2.4 dB. The delay spread 
values for Tx-Rx distances of 2-12 m ranged between 2-3.4 
ns. Furthermore, asymmetric side reflections were investigated 
for scattering components when specular reflections are absent, 
where an additional path loss of 7 dB was observed when the 
vehicle moved by 2 m, and increased up to an additional 10 dB 
for longer distances. Finally, for the multi-lane case of front 
and rear reflection, the reflection losses from vehicles on the 
rear side were measured between 15 to 30 dB, and 24 to 42 
dB from the front side reflection. 

Moreover, Yi et al. [197] examined the Vehicle-to- 
Infrastructure (V2) propagation channel for two urban and six 
highway scenarios at 300 GHz, under clear, rainy and snowy 
conditions based on measurements conducted in [154]. The 
attenuations due to rain and snow in the propagation channel 
were considered to contribute to the path loss by 26 dB/km 
based on the ITU model and 348 dB/km based on measure- 
ment, respectively. For the V2I urban scenarios and for Tx 
antenna heights of 5 and 25 m in LoS with the Rx, mean 
Rician K-factor and delay spread values were recorded as 7.32, 
—0.04 dB and 7.58, 26.62 ns, respectively, where the nega- 
tive correlation between K-factor and delay spread indicates 
the presence of more effective MPCs reflected from build- 
ings on the side of the road at greater Tx antenna heights. 
Furthermore, angular spreads in elevation were found to be 
much smaller than azimuth spreads, indicating that most of 
the MPCs came from the buildings on the side of the road 
which decayed much faster in case of rainy/snowy conditions 
in comparison to clear sky conditions. For the V2I highway 
scenarios were the Tx antenna was varied at 10, 15 and 25 m 
and placed at LoS to the Rx, the effects with and without a 
metallic noise barrier and a 0.8 m brick fence were consid- 
ered. Rician K-factor values of, e.g., 1.71 and 2.68 dB with 
and without the barrier at 10 m show that the reflections have 
a significant influence on the channel. At Tx heights of 25 and 
5 m, values of 4.3 and 2.32 dB with barrier were recorded, 
showing that due to smaller reflection angles at 25 m, there is 
less significant effect in terms of multipath propagation. The 
delay spread, showed the highest values with the presence of 
the barrier at a 10 m height, of 9.11, 8.88 and 6.65 ns for 
clear, rainy and snowy conditions, respectively, demonstrating 
strong reflected components arriving at the Rx, with a signif- 
icant influence of rain and snow droplets. In contrast, for the 
other cases the contributions of rain and snow droplets was 
insignificant. 
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From the above contributions, it can be concluded that parts 
found on a vehicle such as bumpers, and windshields induce 
much greater reflection and scattering losses in comparison 
to the main metallic body of the car, while reflections from 
nearby vehicles were seen to contribute to the overall signal 
strength, especially when found closer to the transmitting vehi- 
cle. Height variations and distance separations of Tx and Rx, 
were seen to affect the signal strength significantly, causing 
delay spread variations and broadened peaks. When the spec- 
ular component from a nearby vehicle is not present, additional 
losses would occur and multipaths would be weakened. The 
effects of snow and rain particles may cause multipaths that 
can contribute to the received power in clear conditions, to 
diminish. Furthermore, the attenuation due to snow does not 
follow the ITU recommendation for rain attenuation, indicat- 
ing that more measurements need to be taken, for a proper 
model to be developed at THz. Furthermore, future works 
need to account for traffic flow causing Doppler shifts, since 
LOS communication cannot be guaranteed between different 
vehicles. 

8) Nano-Network Channels: Jornet and Akyildiz [103], 
presented a capacity evaluation for EM wireless nano networks 
in the THz band between 0.1-10 THz. Three different eval- 
uations were carried out, (i) by varying the power allocation 
scheme, (ii) by varying the molecular composition and (iii) by 
varying the pulse shape. For (i), when the transmitted noise 
Power Spectral Density (PSD) was considered flat over the 
entire band, the scheme performed well at a few tens of mm 
due to negligible absorption losses, tending to the optimal p.s.d 
and reaching capacity values between ~ 1014-1015 bit/s. For a 
case where a transmission window at 350 GHz is used that is 
51 GHz wide, the capacity would be two orders of magnitude 
below the optimal PSD. For (ii), by varying the concentration 
of water vapor molecules for different transmission distances 
for the different allocation schemes as in (i), as well as for a 
100 fs long Gaussian shaped pulse, the results showed with 
an increase in molecules at small distances a small capacity 
degradation, and severe capacity degradation at longer dis- 
tances. For (iii), it was found that for Gaussian shaped pulse 
durations between 50 and 150 fs, that when the pulse width is 
increased, the p.s.d. becomes sharper, with the power mainly 
held within the center of the pulse. While shorter pulses have 
also shown a better capacity performance, this effect helps to 
increase the achievable capacity as a function of increasing 
transmission distance. 

Elayan et al. [80], presented a channel model along with 
a link budget analysis for In vivo Wireless Nano Sensor 
Networks GWNSN’s) in the THz band. Propagation phenom- 
ena inside the human body were analysed mathematically, and 
the path loss was computed while taking into account wave 
scattering from body particles, besides molecular absorption 
from human tissue and the effects of spreading caused by 
the expansion of the propagating wave. The theoretical results 
showed a loss of ~ 80 dB for blood, skin and water samples 
at a 10 mm distance at a frequency of 1 THz. The absorp- 
tion coefficient calculated, showed an approximately linear 
decrease for wavelengths between 0.1 um to 3 mm, with val- 
ues ranging between ~ 200 dB to 70 dB, 80 dB for skin, 
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water, respectively, and ~ 300 dB to 90 dB for blood sam- 
ples. Results for the scattering coefficient show that its effect is 
negligible compared to absorption, as the wavelength is much 
greater than the size of the body particle scatterers at THz 
frequencies, and can therefore be modeled as Rayleigh scatter- 
ing. A link budget analysis was then conducted for intra-body 
communication between two nano-devices considering a trans- 
mit power of 1 mW and a 2.15 dBi gain antenna for a distance 
of 1 mm corresponding to a path loss of 65.8 dB, with an SNR 
value of 10 dB, resulting to a receiver sensitivity requirement 
of —105.8 dBW. 

Yang et al. [203], investigated intra-body propagation losses 
as a function of dielectric loss of human tissues besides 
frequency and distance. It was found that the effect of power 
loss due to layered structures resulting from the interface in 
between, is limited when the permittivity of the layered sam- 
ples is similar, as for the cases of skin and blood samples. 
Furthermore, the effect of noise temperature was also con- 
sidered, and it was found that it increases as a function of 
frequency and distance reaching up to 310 K in the order of 
a few mm, thus, resulting in an increase of the noise power. 
Capacity levels as high as 100 Tbps could be achieved between 
0.5 to 1.5 THz, depending on the propagation environment and 
pulse shape. Furthermore, it can be deduced that for narrower 
Gaussian shaped pulses as well as by increasing the derivative 
order, the capacity increases due to the increase in the p.s.d. 
within the bands under consideration. 

It can be said, that propagation in EM nano-networks is 
dominated by spreading and absorption losses, in comparison 
to scattering losses and permittivity variations in various sam- 
ples which are much smaller. Furthermore, Gaussian shaped 
exciting pulses have been most commonly employed, and the 
performance could be optimally tuned by varying the shape of 
the pulse and pulse duration. Propagation losses have shown an 
almost linear increase with frequency of operation and transmis- 
sion distance, while an increase in concentration of molecules 
has shown to have a major effect on the propagation losses 
experienced in intra-body mediums. Moreover, by considering 
high gain antennas, propagation through the multi-layered body 
structure can be achieved, and the THz wave could penetrate 
different layers such as skin, water and blood. Finally, the THz 
wave may interact with much smaller particles inside the human 
body at higher frequencies avoiding molecular resonances, and 
thus simplifying propagation modeling. 


A. Summary of 6G Channel Modeling Efforts 


Table VIII summarizes and compares the most important 
scenario-specific modeling contributions found in literature. 
Since channel models in the THz band mainly rely on chan- 
nel measurements, many more measurement campaigns must 
be undertaken in order to realize 6G applications and make 
use of the vastly available spectrum between 0.3-3.0 THz for 
future communication systems design. Having a closer look, 
we can deduce that these efforts lack characterization of the 
channel particularly above 300 GHz and are quite sparse for 
some of the envisioned applications. This is mainly due to the 
low attenuation expected in comparison to the higher bands, 


1988 


IEEE COMMUNICATIONS SURVEYS & TUTORIALS, VOL. 24, NO. 4, FOURTH QUARTER 2022 


TABLE VIII 


IMPORTANT MODELING EFFORTS CHARACTERIZING THE TERAHERTZ CHANNEL FOR 6G WIRELESS COMMUNICATION APPLICATIONS 


Ref. Application Scenario Freq. (GHz) Modeling approach Antenna Channel Statistics Measurement 
Campaign 
141 Kiosk download LoS, multiple order reflection 220-340 RT-stochastic SISO Fading amplitude, AoA, AoD, phase, K-factor, VNA 
delay spread 
21] Kiosk download LoS, parallel, tilted, tilted and par- 220-340 RT-stochastic SISO Rician K-factor, RMS delay spread, fading am- VNA 
allel, multiple order reflection plitude, ToA, AoA, AoD 
143 Kiosk download LoS, NLoS, parallel, tilted, tilted 220-340 RT-stochastic SISO Fading amplitude, path loss, BER VNA 
and parallel, multiple order reflec- 
tion 
145 Data center LoS, OLoS 300-312 Stochastic 4x4 MIMO path loss, shadowing gain, RMS delay spread, VNA 
fading amplitude, Rician K-factor 
146 Data center OLoS with stationary and vibrating 300-312 Stochastic 4x4 MIMO path loss, shadowing gain, RMS delay spread, VNA 
cables Doppler shift, Rician K-factor, delay dispersion. 
22] Data center LoS, NLoS, D-NLoS, Rack-to- 300-312 Stochastic 4x4 MMO, path loss, shadowing gain, RMS delay spread, VNA 
Rack, Blade-to-Blade SISO doppler shift, Rician K-factor, delay dispersion, 
coherence bandwidth 
142 Data center LoS, NLoS, specular reflection 300.2-308.2 Stochastic SISO AoA, AoD, ToA, FSPL, path attenuation, polar- sc 
ization 
147 Data center LoS, NLoS, OLoS 300-320 Stochastic SISO Delay spread, path loss, cluster powers, shadow- VNA 
ing gain, correlation 
190 Data center LoS, R-NLoS, OLoS, NLoS 300-320 Stochastic SISO Delay spread, path loss, cluster powers VNA 
144 Data center LoS, NLoS 300-320 Stochastic SISO Delay spread, path loss VNA 
148 Chip-to-Chip LoS, D-NLoS, NLoS 300-314 Stochastic SISO Path loss, delay spread, coherence bandwidth, VNA 
amplitude fading 
23] Chip-to-Chip LoS, D-NLoS, NLoS 300-314 Stochastic SISO Path loss, delay spread, coherence bandwidth, VNA 
amplitude fading 
198 Indoor LoS, NLoS 300 Calibrated-RT 3x3 MIMO amplitude fading, delay spread, angular spread, VNA 
K-factor, vertical, horizontal polarizations 
149 Indoor LoS, NLoS 300 GBSM SISO Interference, amplitude fading, diffraction, ToA, VNA 
AoA, AoD,path loss, FSL, delay spread 
137. Indoor LoS, Specular reflection, scattering 500-750 Stochastic SISO PDP, fading amplitude, ToA VNA 
78] Indoor LoS, NLoS, scattering 350, 650 Calibrated-RT sIsO Clustering, path loss, AoA, AoD, SNR, Capacity, VNA 
delay spread, coherence bandwidth 
154 T2l/inside-station, LoS, Light-NLoS, Deep-NLoS 300-308 quasi-deterministic, SISO path loss, shadow fading, RMS delay spread, sc 
Intra-wagon, I2I calibrated RT Rician K-factor, angular spread, channel gain, 
XPR, cross correlations capacity 
155 T2I/inside-station, LoS, Light-NLoS, Deep-NLoS 300-308 quasi-deterministic, SISO path loss, shadow fading, RMS delay spread, sc 
Intra-wagon, I2I calibrated RT Rician K-factor, angular spread, channel gain, 
XPR, cross correlations capacity 
156 T2I/inside-station, LoS, Light-NLoS, Deep-NLoS 300-308 quasi-deterministic, SISO path loss, shadow fading, RMS delay spread, sc 
Intra-wagon, I2I calibrated RT Rician K-factor, angular spread, channel gain, 
XPR, cross correlations capacity 
157. T2I/inside-station, LoS, Light-NLoS, Deep-NLoS 300-308 quasi-deterministic, SISO path loss, shadow fading, RMS delay spread, sc 
Intra-wagon, I2I calibrated RT Rician K-factor, angular spread, channel gain, 
XPR, cross correlations capacity 
133 V2V LoS, reflection, scattering 300 deterministic, SISO Multipath interference, SNR, link capacity, SIR, THZ-TDS 
stochastic reflection and scattering losses 
158 v2v LoS, side reflection, rear reflection, 300 NGBSM 4x4 MIMO, Path loss, path gain, AoA, AoD, ToA, delay Ne 
under-vehicle propagation, block- SISO spread, angular spread 
age 
197 V21 LoS, NLoS 300-308 Calibrated-RT SISO Path loss, Rician K-factor, AoA, AoD delay sc 
spread, angular spread, rain, snow 
103 Nano-network LoS 0-10000 Stochastic SISO Path loss, molecular absorption noise, capacity, 
power allocation, radiative transfer theory 
79] iWNSN LoS, scattering 100-10000 Stochastic SISO Intra-body Spreading, molecular absorption, scat- 
tering loss by human tissue 
203 Intra-body LoS 500-1500 Deterministic, SISO Capacity, path loss related to frequency, distance 
stochastic and dielectric loss of human tissue 


as well as because of equipment availability, since the cost 
to purchase such equipment is quite high and also due to the 
limitations of the equipment available to measure much above 
300 GHz. The characterisation efforts performed so far that 
revolve around 300 GHz as well as current contributions for 
the spectrum above, can therefore be used as a basis provid- 
ing a useful comparison for future channel characterization 
in the THz regime. Furthermore, it is obvious that there is 
a lack of sounding equipment with MIMO capabilities that 
are imperative for realizing long distance communications and 
achieving high throughput levels. A more thorough discus- 
sion regarding challenges and potential solutions as well as 
research directions for future channel modeling efforts follows 
in Section VIII. 


VIII. CHALLENGES AND FUTURE RESEARCH DIRECTIONS 


From the work analyzed in this survey paper, there are some 
open issues and challenges identified in measurement and 


modeling of the propagation channel that need to be addressed, 
paving the way for potential future research directions and 
practical solutions to be explored. 

1) Diffuse Particle Scattering: The issue of diffuse particle 
scattering at THz frequencies, where we get rough surfaces due 
to the Rayleigh criterion, becomes crucially important when it 
comes to accurate modeling for use cases applicable in indoor 
environments and within devices. Multipath resulting from the 
layers in the ‘smooth looking’ composite materials as well 
as interactions where the surface height varies the order of 
the wavelength, cause cluster formation and broadened peaks, 
which are some of the effects caused when surfaces with dif- 
ferent roughness levels and dielectric constants interact with 
the THz wave [96]. These effects, cause variations in the 
received power levels from point to point and over differ- 
ent angles, as well as increased delay spreads that may affect 
system design as observed from the analysis conducted in 
Section VII, and it is therefore of great importance to study 
these effects for accurate channel modeling. 
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2) Multipath in the Illuminated Region: The beamwidth of 
the highly directional horn antennas employed and what mul- 
tipaths it picks up in the area it is illuminating, are key factors 
that have not yet been explored when it comes to multipath 
scattering. These factors may cause effects such as unpre- 
dictable multipath arrival and varying frequency selectivity. 
Therefore, evaluating the relation between multipath scattering 
and beam spot diameter at various points and angles, becomes 
increasingly important for channel modeling at THz, requiring 
an in-depth investigation. 

3) Capturing Time-Varying Channel Characteristics: So 
far channel characterization mainly focused on static scenarios 
where PtP links will be employed. But in reality most cases 
will have to adapt to dynamic variations of the propagation 
channel due to the highly directive characteristic of the THz 
waves. Understanding the characteristics of the reflections off 
surfaces becomes important, whether it comes from a human 
obstructing the LoS path or if data is received from a vehicle 
that is passing by and there is a NLoS bounce that varies as 
the vehicle is moving. This differs from the traditional Doppler 
spread modeling, and modeling of these effects could benefit 
by an AI approach that can relieve the need for RT updates 
based on previous data. Similar to the approach employed 
in [204] that allows for dynamic AoA estimation to cater for 
human movement in THz links by adjusting an adaptive direc- 
tional antenna. By optimizing the state transition probabilities 
through knowledge of previous estimates, this AI approach has 
observed to outperform RT and human movement modeling. 
In a similar manner, approaches as such, can be employed to 
model scatterer mobility, as well as surfaces found in the envi- 
ronment from point to point over wider bandwidths, extending 
limited data sets, and solving the problem of material library 
formation employed in commonly used RT tools where the 
material properties are considered to be constant, altogether 
with manageable complexity. 

4) RF EMF Exposure: The expected propagation ranges 
specified in Table III, indicate that many applications scenar- 
ios are confined to small propagating environments, where the 
user may be subject to long term exposure of THz radia- 
tion. Such scenarios include Intra-wagon links and wireless 
connectivity in indoor locations. These applications could 
potentially offer Virtual Reality (VR) services, where expo- 
sure effects on the body may be different, more in regard 
to where the devices may be used on the body, such as in 
VR headsets that then could directly couple or radiate to the 
skin and eyes. Furthermore, beams of high EIRP are required 
for links in enclosed areas, and it could be frequently the 
case that a person moves into the beam for a short or long 
period of time. The International Commission on Non-Ionizing 
Radiation Protection (ICNIRP) has so far defined exposure 
guidelines up to 300 GHz [205]. Subsequently, it becomes 
a case of what happens beyond that range that will allow 
limits to be set. THz radiation is generally considered as a 
form of non-ionizing radiation [206]. Recent studies however, 
have shown that with high intensity transmission, different 
biological molecules exhibit absorption spectra, inducing res- 
onances that may damage but also repair DNA and protein 
cells [207]. Research is still at its early stages and due to 
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the lack of measurement data, we are yet unable to determine 
THz health effects to establish radiation limits caused by near 
field short-term and more importantly long-term exposure, to 
enable guidelines to be written indicating exposure limits for 
safe use of THz radiation. 

5) Enhanced Channel Simulators: Current simulators, fail 
to address the problem of accuracy and complexity, due to the 
need for extensive measurements and long simulation times 
to obtain high accuracy in the results obtained. Furthermore, 
these simulators fail to consider dynamic channel variations 
that result from interactions with scatterers found in the prop- 
agation environment. The channel simulators developed, must 
also be able to cater for synergy with the lower bands, where 
REF/THz and THz only operation can be expected, such as 
for fronthauling, and allow for dynamic spectrum allocation 
schemes, e.g., for multiple access techniques operating under 
the NOMA principle, which are expected to prevail in 6G 
coexisting networks. This however would induce a great deal 
of additional complexity due to the exhaustive computation 
of the Beer-Lambert law, that has to be managed efficiently. 
Furthermore, current channel simulators do not consider oper- 
ation with RISs which is vital in improving the coverage in 
6G networks, and therefore must be integrated. 

6) RIS Assisted Channel Modeling: RISs, hold great poten- 
tial for application at THz to direct wave propagation when 
the LoS path is not available. Therefore, RIS assisted chan- 
nel modeling must be considered in PDP type models. SV 
type modeling can be incorporated for cluster identification, 
where the clusters formed are expected to be much stronger 
in comparison to reflections without the RIS, and then it is a 
case of how the MPCs within the clusters are modeled, and 
whether commonly employed mixture models can also apply 
to describe them effectively. This would benefit from larger 
bandwidths offering higher delay resolutions, and the ability 
to distinguish multipaths that are densely packed around the 
peak of the main delay component. Therefore, modeling RIS 
propagation in comparison to reflections without RISs, and 
how to incorporate existing clustering models as well as AI 
approaches to capture channel non-stationarities, gives way for 
more work to be done toward that direction. 

7) Hybrid Channel Modeling: It is apparent that the hybrid 
modeling approach where empirical data are fed to RT simu- 
lators to extend the limited data sets is very popular. Similarly, 
channel characterization should be expanded and comple- 
mented by the use of various stochastic and deterministic 
models to describe each part of signal propagation from one 
end of the communications system to the other. As an exam- 
ple, many empirical data sets can be collected and an SV 
type model can be used for grouping MPCs together and a 
second order GMM can be used to describe MPCs within 
arriving clusters. Furthermore, FDTD/MOM modeling can be 
employed to model a portion of the propagating environment 
when RT simulations do not provide enough information about 
the channel. It provides much higher accuracy and manages 
computational complexity since it can be applied to model 
areas of complex discontinuities [176]. These methods can be 
then tuned within a unified model in parallel with great knowl- 
edge of the weather conditions and scatterer properties, and 
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employed along with ML techniques that can capture dynamic 
channel variations, and which could perhaps allow for this 
characterization to be expanded over larger bandwidths, result- 
ing in reduced computational resource requirements. This is 
due to the difficulty and timeliness of collecting a large amount 
of useful data using a VNA, and due to technological limita- 
tions of real-time channel sounders which are preferred for 
collecting data, but are limited by their short bandwidths and 
frequency scaleability. 

8) Enhanced Channel Sounding Systems: The distinct 
propagation features of the THz wave as well as the unique 
propagation environments expected for the envisioned appli- 
cations, give rise to distinct measurement requirements which 
have been explained in Section III. The very high bandwidth 
availability, requires channel sounding systems to operate over 
a large bandwidth and at a scalable center frequency, some- 
thing that is currently available with VNAs and spectroscopy 
systems. The high carrier frequencies, also indicate that anten- 
nas operating at these frequencies will result in highly directive 
pencil beams. This brings forward the issue of channel spar- 
sity, that would require a high dynamic range of the system 
in order to detect strong enough signal components result- 
ing from LoS as well as multipath propagation due to the 
large losses present in free space and the atmosphere, that 
would enhance the overall quality of the signal received. It is 
therefore apparent that channel sounding systems with MIMO 
capabilities are in demand. Furthermore, the formation of 
pencil beams, coveys that in order to efficiently character- 
ize propagation with limited amount of information, it would 
be required to have a very high multipath resolution, thus 
requiring a large bandwidth such as that offered by VNAs 
and spectroscopy systems. Furthermore, the channel sounding 
systems developed, should ideally operate under the SC princi- 
ple, that measures the samples without surpassing the channel 
coherence time, to allow Doppler shift variations to be cap- 
tured to analyse dynamic effects such as link blockage. It is 
also evident, that channel characterization over long distances 
is currently limited by phase synchronization of the signal. 
Therefore, when sending a signal from Tx to Rx, we can only 
measure delay paths relative to the strongest component of 
the PDP. Thus, the need for reference signals to maintain syn- 
chronisation over short distances with the use of lossy RF 
cabling, or otherwise, high precision rubidium/caesium clocks 
on either side to replace the synchronization cable, which are 
generally accepted practices. However, by adopting the much 
more efficient RFoF technique, measurements could be con- 
ducted at much greater ranges with negligible synchronisation 
errors which comply with timing standards. Therefore, high 
performance channel sounding systems are needed, that can 
maximise the above mentioned capabilities, to eventually help 
improve channel characterization efforts. 

9) Extensive Measurement Campaigns: From the review 
carried out in this work, it is evident that while certain use 
cases and scenarios have been characterized sufficiently in 
the 300 GHz band, some other ones lack important char- 
acterization efforts. Moreover, there is a great need for 
further characterization works to be conducted within the 
atmospheric windows identified in Section IV, before further 
standardisation processes take place. Furthermore, new use 
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cases and application scenarios that are appealing in the THz 
bands and have yet to be characterized, such as communi- 
cations in ships, airports, and Device-to-Device (D2D), also 
require extensive measurements. Nonetheless, other applica- 
tions such as Unmanned Aerial Vehicle (UAV) communica- 
tions and Inter-Satellite Links (ISLs) may not be measurable 
as the system is not terrestrial and channel characterisation 
may mainly rely on simulations. 

10) Empirical Data Acquisition: As empirical path loss 
models do not exist, and modeling of the propagation channel 
heavily relies on interactions with the propagation environment 
that cause unpredictable scattering and frequency selectivity, 
large data acquisition is essential. Data collection using high- 
precision mechanical rotational and translational units that are 
required in order to achieve a high autocorrelation is time 
exhaustive, where, e.g., data collection for angular-dependent 
PDPs for different positions, may take 100s of hours [194]. 
A practical solution to this problem, which nonetheless has 
been underutilized at THz due to the high cost of applica- 
ble positioning equipment and complexity in setting up its 
operation, is to employ General Purpose Interface Bus (GPIB) 
control [208]. GPIB control, would allow the antennas to be 
varied in the spatial and/or angular domains, and the data 
from the VNA to be saved automatically, thus, enabling an 
extremely large set of data to be acquired in a much shorter 
time period. This option seems the most beneficial solution 
for large empirical data acquisition given the current status 
of dynamic channel sounding technologies that are limited to 
operation at 300 GHz, and the only account of this system in 
literature for channel measurements at THz, is reported in [79]. 


IX. CONCLUSION 


In this paper, we provided an overview of propagation 
phenomena, measurements and modeling of the propagation 
channel at THz. An introduction to the THz band was ini- 
tially given, discussing enabling technologies, standards and 
potential applications that could be supported, while exam- 
ples of measurement requirements for realizing the envisioned 
applications were provided. THz propagation phenomena such 
as path loss and molecular absorption as well as propaga- 
tion mechanisms were then introduced, and possible bands 
of operation were assessed for 0.3-3 THz. THz channel 
sounding methods were described and assessed in relation 
to measurement requirements for the applications identified, 
and existing campaigns found in literature were presented. 
Modeling approaches with potential applicability at THz for 
characterizing propagation for the envisioned applications and 
current modeling efforts were evaluated. The unique charac- 
teristics of the propagation channel for each of the envisioned 
applications were then presented and analyzed. Finally, an 
insight was given to the issues and challenges faced for mea- 
surement and modeling in the THz spectrum, and potential 
solutions and research opportunities that must be explored to 
tackle these challenges were discussed. This paper provided 
a solid understanding of channel propagation at THz, setting 
a good knowledge base for future work to be conducted to 
realize various applications targeted at this band for operation 
in 6G networks. 
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